On Classical Electromagnetic Fields 



i. preliminaries: a review of some basic concepts and methods: 
The Microscopic Maxwell's Equations in the Time Domain : 

VxE(r,t)=-yB(r,t) [I-la] 
VxB(r,?)=|io j(r,?)+8o|io|^E(r,r) [Mb] 



V ■E(r,t) = -p(r,t) [I-lc] 



VB(?,r) [Md] 

where the meaning of fields E(r,?) and B{r,t) is ultimately defined in terms of the 
Lorentz force on a charge q moving at a velocity u — viz., 

F(r,t) = qE(r,t) + quxB(r,t) [1-2] 

The Macroscopic Maxwell's Equations in the Time Domain 

Following common practice, we set forth a particular set of macroscopic Maxwell's 
equations that applies in the high frequency or optical regime.^ 

VxE(i,t)=-yB (?, t)=-iioj-H (i,t) [ I-3a ] 

V xB (r,?) =|i V xH (r, t) = [Iq (f , + ^ P (r, r) j +£ ^ [Xq ^ E (r, t) [ I-3b ] 



1 All bound current effects are included in the polarization density, since magnetization density "ceases to have 
any physical meaning at relatively low frequencies." See Section 62 in L. D. Landau and E. M. Lifshitz, 
Electrodynamics of Continuous Media, Pergamon Press (1960). 
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V ■£(?,?) = - [p(?,r)-V ■?(?,?)] [I-3c] 
V.B(f,?)=^„ V H (?,?)= [i-3d] 

A Phenomenoloical Representation of the Linear Dielectric 
Response of Matter ;^ 

The following is the most general phenomenological representation of the linear 
dielectric response of a given material that incorporates dissipative, non-local, and 
anisotropic effects: 

In a tensor representation 

(r,t) =8o j dt'\\\ t/) Ep (r, t) dv [ I-4a ] 

In a dyadic representation 

t 

P'"' (r, t) = j J^jj % (?' f ^' -E (r, df' [ I-4b ] 

For the present and for most of our discussions, we neglect nonlocal effects and treat 
only dispersive (dissipative) and anisotropic effects so that 



2 In our later treatment of nonlinear optics, we will begin by adding the following nonlinear phenomenological contributions: 
' Eo E J J J J df,dt,df,dt, t- t,;f- ?2. t- t,) Ep(?, E^(?2, t,) 



+£oILJJJJJJ df^dt^df^dt^dfjdt^ Xa?,js t- t{, r- r^, t- t{,T- t- t 



3/ 



PT^ I, i2 I3 /J /2 '3 
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P{r,t)=£, \ dfx{f,t-f) E{r,f) 



[1-5] 



or p(r,a))=8o% (r,co) E(r,o)) [1-6] 

t oo 

where % (r, co) = j % (r, f) exp [-z (£>{t- t')] df = ^ x (? . ^) exp [-/(OX ] [ 1-7 ] 



IMacroscopic IMaxweirs Equations in the Frequency Domain Valid for Linear. 
Local, Anisotropic IMedia in the Optical Resime. 



V X E (r, (o) = V X E (r, co) • D (r, co) = -id) B (r, (o) = IJ-q H (r, O)) 

V X B (r, 0)) = |io V xH (f , (o) = [ig J (f, (o) -I- i(0 [i^ 8 (f , (o) • E (r, cd) 
= IJ-o J (r, (o) -I- /(o |io D (r, co) 

V e (r, (o)- E (r, co) = V D (f , co) = p (f , co) 
V ■B(r,co) =1^0 ^ ■H(r,co)= 



[I-8a] 
[I-8b] 

[I-8c] 
[I-8d] 
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Helmholz equations for the Frequency Domain Vector and 
Scalar Potentials in a Uniform, Linear, Isotropic Dielectric 

We define the (Magnetic) Vector Potential as 

H(r,0)) = — B(r,co)= — VxA (r,0)) [1-9] 

which, by design, automatically satisfies one of Maxwell's equations — viz., [ I-8d . 
We introduce the (Electric) Scalar Potential in the form 

E(f,co)=-/co A(f,a))- V^(r,co) [I-IO] 

which, again, automatically satisfies another Maxwell equation — viz., [ I-8a ] 
Therefore, for uniform, isotropic media Equation [ I-8b ] becomes ^ 

V xV xA (f,(o) =\Xq j(r,a))-i-(oV,)£(a)) A (r,(o)- zo)|io£((o) V(|)(r,(o) 

rills.! 

= VV-A(r,(o)-V'A(r,(o) ^ 

and Equation [ TSc ] becomes 

-i(oe{(o) V A (r,(o)-£ (co) V^(|) (r,(o)= p(r,(o) . [ i-Hb ] 

Since 

V A (r,co) is as yet undefined, we define it in the Lorentz sause as 

V A (r,a)) = -/a)|j.oe((o)(^(r,a)) [1-12] 



3 Since div curl { }=VVx{ }=0. 

^ Since curl grad { } = VxV{ } = 0. 

5 Using ax (bxc)= b(a c)-c(a- b). 
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to simplify Equations [ I-l la ] and [ I-llb ]. Thus 

V' A (r,o))+0)'|ioeN A (f,(o)=-|io J (r,a)) [I-13a] 

and 

V^(^ (r,(o)+a) Vo <^ (f ,co) = (r,(o) [ I-13b ] 

Therefore, in the Lorentz gauge both A (f ,co) and (? ,£o) satisfy inhomogeneous 
(and homogeneous) Helmholz equations! 

However, in the Coulomb muse we define V ■ A(r,a) ) = and then Equation [ I-8c ] 
becomes 

V>(?,a))=-^p(?,(o) [I- 14a] 

Conservation of charge requires that V J (f ,co)- /cop (r,co) = so that Equation [ I-l la ] 
becomes 

V' A (r,(o)+coVoe N A (r,co) =-^1^ J^'™^' (?,a)) [ I-14b ] 

where 

J<'"-^ M=j(?,co)-j('°"^'(f,a)) 



J (r,co)- [-/C08 (co)V(|) (f ,0))] 



[1-15] 
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II. RAYS: THE EIKONAL TREATMENT OF GEOMETRIC OPTICS 6 

Since ancient times, the notion of ray or beam propagation has been one of the most 
enduring and fundamental concepts in optical physics. As a zeroth order approximation we 
might consider a plane wave to be a model of a beam and its propagation vector to be a 
model of a ray. This is a reasonable start, but it is a much too restricted view and we can 
do much better. What we need is a solution to Maxwell's equations which is like a plane 
wave, but limited in spatial extent. One approach, the simplest, is called variously ray, 
Gaussian or geometric optics. 

A Ma xwellian Derivation of the Eikonal Equation: 



To fully understand geometric optics in the context of Maxwell's equations 
writing the electric and magnetic fields as pseudo-simple waves — viz. 



we start by 




[Il-la] 




[Il-lb] 



where ko= ©.^ii^ =(0 /c 



It is assumed that 

e(f,(o) and h(r,co) 

are weak functions of position . 

phase function S(r,(o) is the spatially varying phase of the pseudo-simple 
the cases of pseudo-plane waves and pseudo-spherical waves the phase 
given, respectively, by 



wave. For 



function is 



The scalar 



koS(r,(o) = A:k^-l-3;k^,-l-zk, 



[ n-2a ] 



6 



See, for example. Max Born and Emil Wolf, Principle of Optics, Pergamon Press (1986), Chapter 3. 
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and koS(r,a)) = + 7 [ II-2b ] 

We now substitute these pseudo-simple wave expressions (i.e. Equations [ II-l]) into 
Maxwell's equations to obtain 

exp[-/koS (?,«)] I V xe (f,co)- iko VS(f,a))xe (f,a))|=-/|i.oCk(,h (f,a)) exp[-ikoS(f,a))] 

exp[-jkoS(f ,co)]|vxh(f,(o)- ikg VS(f,(o)xh(f,(o)|= i8(f ,(o)ckQe(f,co) exp[-jkoS(f ,(o 

exp[-jkoS(f ,(o)]|v |E(f,(o)e(f,(o)]- jk(,E(f,(o) VS(f,(o) e(f,(o) j= 

exp[- jkoS(f ,0))] {v •h(f,(o) - ik^ V S(f ,(o) • h(f ,«) } = 

Rearranging, we obtain 

VS(r,a))xe(r,a))-|a,oC h(r,a))= [/kg] ^Vxe(r,a)) [n-4a] 

VS(f ,0)) X h(f ,0)) +8(f,a)) ce(f ,0)) = [ik^]^ V xh(f,a)) [ II-4b ] 

VS(r,to) (r,a))e(r,a))] = [ik^y' V [8(r,a))e(r,to)] [ n-4c ] 

V S (r,a)) • h (r,a)) = [ i k^ ] V h (r,©) [ II-4d ] 

In the ray, Gaussian or geometric approximation we assume that we may neglect 
the RHS's of these equations. To get something useful we multiply through the first 
equation (i.e. Equation [ II-4a ]) as follows: 

[|ioc]"'vs(f,co)x{ Equation [ II-4a ] } [n-5a] 
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[iIqc] ^VS(r,a))x{vs(r,a))xe(r,a))-|ioCh(r,(o)}«0 [II-5b] 

Applying the "abc = bac - cab" rule^ we obtain 
[^lo c]"'{ V S(f,a)) [ V S(f,a)) • i(f,to)] 

-e(f,(o) |vs(r,co)|'-|iocVS(f,a))xh(f,(o)[ ==0 



[n-5c] 



which becomes upon substitution from the second Equation [ II-4b ] 
[|ioc]"'{vs(f,(o)[vs(f,(o)e(f,(o)] 

- ,0)) |v S(r,o))| i+8(r,o)) ce(f ,co) ~ 



[ II-5d ] 



From Equation [ II-4c ], we see that the first term vanishes in the geometric 
approximation — i.e., if we neglect the term [^k^ ] V je (r,co) e(r,(o)] . Therefore, 
for non- vanishing e(f,03) we obtain the following reduction of Maxwell's 
equations : 

VS(f,(jo)|' ==E(f,(o)|io c' = n'(f,co) [ II-6 ] 

where '^(f ,co) is the index of refraction. More explicitly, we may write an equation for 

a "ray vector" ~ i.e. the tangent to a space curve orthogonal to the surfaces of constant 
S(r,(o) 



Again using ax(b X c) = b(a • c) - c(a • b) . 
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VS(r,a)) = n(f ,0)) s = ^(r,©) — 

ds 



We illustrate the geometric relationships below: 



[11-7] 





Reference poin't 



grad S(r^w) 



We may now derive the all important eikonal equation. To that end, we first take a derivative 
along the ray direction - viz. 



— [vs(f,co)l= — n(r,(o) — 



ds 



ds 



[ II-8a ] 



However, from the definition of the grad operator we know that 

i/[vs(?,a))] = dY-V [vs(?,a))] 



so that 



d 



df 



f • vfvs(?,co)l = 1-^ VS(f,a))|. Vfvs(f ,co)l = -f n(f ,«) ^ 
ds ^ ^ [«{r,toj J ^ dsl ds 



[ II-8b ] 
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or 




1 



v{vs(f,a))VS(f,a))} 




1 



v{n'(r,a))} = Vn(r,a)) 



J r \ dr 
= — n r,a) — 
dsl ^ ^ ds 



[ n-8c ] 



Thus we have obtain the eikonal ^ equation for the ray vector ~ viz. 



d \ dr 

— n[r,G)) — = 

ds\_ ds _ 




[11-9] 



First Application of the Eikonal Equation: Mirages 

Air adjacent to a hot surface rises in temperature and becomes less dense. Thus over a 
flat hot surface, such as a desert expanse or a sun drenched roadway, air density 
locally increases with height and the average refractive index may be approximated 
by a simple linear variation of the form 



where x is the vertical height above the planar surface, ng is the refractive index at 
ground level, and K is a positive constant. 

We may use the eikonal equation to find an equation for the approximate ray 
trajectory — i.e. an equation for ray height x as a function of ground distance z -- of a 
hght ray launched from a height Xq and at an angle with respect to the surface of the 

earth. 



n {x)= {1+ K^} 



[11-10] 



The eikonal (from the Greek: eiKOSV means image) was introduced in 1895 by H. Brans. 
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Therefore, 



d_ 
ds 



/- ^ dr 
n r,0) — 

' ds. 



= V/i(f,a)) 



d^x 1 



dz n[x,Gi) dx 



d , . 
— nyxisy) 



[n-lla] 



or from Equation [ II- 10 ] 



d^x 
Hz' 



[Il-Ub] 



Thus, the ray trajectory is given by 



f(z) = j^-^ z^ + tanSg z+Xqj x+z 



[11-12] 



Ray trajectories diverted by a hot surface 




20 40 60 SO 100 

Distance along hot pUnar surface 
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Second Application of the Eikonal Equation: The "ABCD" Ray 
Matrices - A Systems Approach to Optics 



1. Uniform dielectric medium — i.e. Mj,(xi) is a constant so that — 

as 



dr 

. ds . 



0. 



Thus, the ray must be a straight line which may be written r = 5 a + b . 

In the two-dimensional paraxial approximation, we assume that s~ z and write 



J dp 

Pout =Pm + L — 

dz 



and ^ 

dz 



dp 

dz 



[11-13] 



where p = p0=xx + yy 



Pin 




^P'out 

Pout 


_- 


^ in 

^in 



1 


— 


^ out 



We may write results of this sort in the form of the famous and highly useful ray 
transform or ABCD matrix -- viz. 



Pout 




A 


B 




Pout. 




c 


D 





[11-14] 
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In the case of a uniforai dielectric 



Pout 




1 


L 




Pout. 







1 





[11-15] 



so that A = 1, B = L, C = 0, and D = 1 



2. A dielectric discontinuity : Starting with Equation [ II-7 ] and noting, once again, 
that curl grad { }=VxV{ } = Owesee that 

VxVS(f,a))= Vx|/i(f,a)) s}= [11-16] 

which is identical to the saltus condition on the electric and magnetic fields at a 
dielectric interface! Hence | n\r,(iij s [ continuous across the dielectric 

J tangent 

boundary so that sin6j = sinGj — i.e. Snell's law! This result in the paraxial 
approximation ( i.e., sinG » tanG = G ) may be written in ray matrix form as 



Pout" 




"A 


b' 




"pin" 




"1 




"P,n" 


Pout. 




c 


D 




.pL 




.0 "in/ "out. 




.pL. 



[11-17] 
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3. A " Thin" lenses : In passing we note that the ray matrix of a thin lens is given by 
or, perhaps more accurately, a thin lens is essentially defined by 



[11-18] 



Pout" 




'A 


B' 




"P,n" 




1 




Pin' 






c 


D 















4. Axially symmetric GRIN media : Consider the use of GRaded INdex technology to 
obtain an axially symmetric variation in the index of refraction of the form ^ 



n{p)= 



r 



l-A 



rp 



'1 



[11-19] 



A note on GRIN technology : In GRIN technology one builds up a glass rod with a specific radial index of 
refraction distribution by fusing a sequence of coaxially arranged glass tubes with appropriate index and diameter 
as illustrated in the following: 

Coaxial dielectric (glass) tubing 
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Within such a GRIN rod, we write f = pp + z z for the ray coordinates and 

V«(f) = (3 — n(p) for the index variation. Using the eikonal equation — i.e. 
dp 

Equation [ II-9 ] ~ in the paraxial approximation, we find 



dz 



[ II-20a ] 



or 



d . s J r . . dp~\ 

nip) ~ — I n p — I 

dp dz[_ ^ ' dz\ 



[ II-20b ] 



Therefore 



— 2 P ' / \ 
dz n\p) dp 



^ :^"(p)=-^i4«(p)] 



[ n-20c ] 



or 



:p- 



dz^ n{p) dp y 



1- A 



fpTll 1 r mAfpf-'l 
\a J \ n\p) a \a J 



niA 



a 



/ \m— 1 
P 



[ II-20d ] 



Doubtless, the simplest and most valuable instance is m = 2 ~ i.e. what is usually 
called parabolic or quadratic material ~ wherein 



d^ 2A 

— P- 

a 



dz 



^pf 2A 

- =- — P=-K P 

Va y a 



so that 



P (z) = Pin cos(kz) + ^ sin(Kz) 



[11-21] 
[ II-22a ] 
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P\z) =-PinK sin(Kz)+p; cos(kz) 
In tenns of a ray transform matrix 



Pout 




A 


B 




J^out. 




c 


D 





r cos(kz) k ^sin(Kz)1 
[-Ksin(Kz) cos(kz) J 



where K= yjlAfi 



Ray trajectories confined in a GRIN rod. 



[ II-22b ] 



[ 11-23 ] 



Imaging Distance 
K AZ = TL 




Focusing Distance 
K AZ = n/2 
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Alternative (Hamiltonian) Derivation of Eikonal Equation: 

FERMAr S PRINCIPLE 

Like most laws of physics, the equations of geometric optics can be derived from a 

variation principle. In this context the variation principle is called the Fermat 
principle which states that a ray always chooses a trajectory that minimizes the optical 
path length - viz. 

J n {x,y, z) ds = minimum [ 11-24 ] 

Pi 



See, for example, Dietrich Marcuse, Light Transmission Optics, Van Nostrand Reinhold (1972). 

11 More precisely, the path must be a local extremum and in rare cases may, in fact, be a maximum. See R. Y. 
Luneberg, Mathematical Theory of Optics, University of California Press, Berkeley and Los Angeles (1964). 
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where the line element, ds , is measured along a ray and the two end-points Pi and P2 
are fixed in space. 12 Analysis of the variation problem is facilitated by choosing the 
projected coordinate z as the new variable of integration. Accordingly, 

ds= -Jdx^ + dy^ + dz ' = yjl + x^ + y^ dz , [ 11-25 ] 

, dx , dy 

where x = — and y = —, Fermat's variation principle is transformed into the more 
dz dz 

familiar Lagrangian form ~ viz. 

Pi 

I L{x, y, x', y')dz = minimum [ 11-26 ] 

Pi 

where L{x,y, x',y') = r^x,y,z) ■^l^^x^^'^^Y^ . [11-27] 

The minimization procedure is then well-known in the variational calculus and 
leads to the famous Euler-Lagrangian equations ~ i.e. 

d dL dL „ 

— —,- — =0 [n-28a] 

dz ox ox 

d dL dL ^ 

— —,-—=0 [II-28b] 
dz oy ay 

When applied to the Fermat Lagrangian, as defined in Equation [ 11-27 ], these 
equations yield 



1^ From Equation [ 11-7 ] we see that 

S{F,)- S{F,)= J n{x,y,z)ds . 
Pi 
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= ^\ + x^ + y^— [II-29a] 



ny' r Ti^n 

=d===^l + x +y — . [II-29b] 



dz ^i+x'^+y^ ^ dy 

Using Equation [ 11-25 ] we see that the Euler-Lagrangian equations may be 
expressed in the vector form as 

[jT dx 1 d I dy 1] _ \ d n d n] nSO 

which is precisely the content of Equation [ II-9 ] -- QED. 

Hamiltonian Formulation of Ray Optics 

The analogy between ray optics and particle mechanics is most striking when the 
equations of ray optics are expressed in Hamiltonian form.i^ To that end, we define the 
generalized momentum which is canonically conjugate to x and y by the vector 
equation 

The Hamiltonian is then define in terms of the generalized momentum by the relation 

H(x,y,p,,Py)= p,x+Pyy-L{x,y,x,y) . [11-32] 

With the assumed functional dependence of the Hamiltonian, we form the derivatives 



1^ The formal theory of optical systems was developed by Sir W. R. Hamilton in 1828-37. 
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dH , dx' dy' dL dx' dL dy' 

^=^+P.^ + Py^ T~^T~ [II-33a] 

op, dp^ dp^ dx dp^ dy dp^ 

dH dx' , dy' dL dx' dL dy' ^ „ , 

^=P^^+y+Py^-^'^ T-^T^ • [II-33b] 

dPy dpy dpy dx dpy dy dp^ 

Given the definitional relationships embodied in Equation [ 11-31 ] we see that these 
expression reduce to one set of Hamilton's equation ~ viz. 



\dx_ dy\_ \dH_ dH [ 
[dz'dzj [dp/ dp J' 

The other set of Hamilton's equation ~ viz. 

\ dp^ ^1 

{dz'dz] [ 3x ' dy 



[11-34] 



[ 11-35 ] 



follow directly from the Euler-Lagrangian equations ~ i.e. Equations [ II-28a ] and 
[ II-28b ] ~ and the definitions embodied in Equation [ 11-31 ]. Using the Fermat 
Lagrangian we see that 

r 1 f 3L 3L 1 [ nx' ny' 1 

{p.,Pyh\ = / . , [11-36] 

[dx dy ] yj\+x'^+y^ + + j 
and consequently that we may solve for { x',y' } in terms of { p,,Py } as 

^ P P ^ 

{-VIH ■ ^ , \ [11-37] 

[ -P.-Py V" -P.-Py J 
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Substituting into Equation [ 11-32 ], we find an expression for the Fermat or ray 
optics Hamiltonian ~ viz. 

which resembles the mechanical Hamilton of a relativistic particle ~ i.e., 



But the analogy is even stronger in the paraxial approximation where the Hamiltonian is 
approximated by an expression which is identical in form with the Hamiltonian of a 
non-relativistic particle ~ viz. 

rr I Px+Py Px+Py r XT 



when and « {n) 



Applying the quantization rules of quantum mechanics to these Hamiltonians, we can go full circle and recover 
wave optics from ray optics. Equation [ 11-38 ] leads directly to the equivalent of the relativistic Klein-Gordon 
equation while the equivalent of the nonrelativistic Schrodinger equation follows directly from Equation [ 11-39 |. 
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III. THE PARAXIAL WAVE EQUATION -- PROPAGATION OF GAUSSIAN 
BEAMS IN UNIFORM MEDIA 
Derivation of Paraxial Wave Equation; 

In point-to-point communication, we may think of the electromagnetic field as propagating 
in a kind of "searchlight" mode ~ i.e. a beam of finite width that propagates in some 
particular direction. In analyzing this mode of wave propagation, we make use of an 
important solution to the so call paraxial approximation of the electromagnetic wave 
equation (or, more precisely, the paraxial approximation of the Heknholz equation). 

To that end, we first derive the paraxial approximation and then examine the free-space 
Gaussian Beam solution(s). We start with the homogeneous Helmholz equation for 
the vector potential in the form — see Equation [ I-13a ] 



We are looking for a wave propagating in, say, the z-direction, so we write a particular 
component of the potential in the form 



The function ^(r,a)) represents a spatial modulation or "masking" of a plane wave 
propagating in the z-direction. The z-direction is obviously special and it is useful to 
appropriately parse the differential operators. For the grad operator we may write 



V' A(r,to)-i-to'|io8(a)) A(r,to)=V' A(r,a))-i-A:^ A(f,to) =0 



[HI-l] 



A„ (r,a) ) =^ (r,o)) exp(-/^z) 



[ni-2] 



grad { 



}=V{ 



}=V.{ 




} [ni-3] 



where, for example. 
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v,{ 



}=xf { 

O X 



By 



so that 



[111-4] 



V A„ (f ,0) ) = V, (r,co ) + zy'^ (^'^ )-ikz r,a) ) j exp{-ik z) [ ni-5 ] 



For the Laplacian operator we may write 



A„ (r,a)) =V ' ^(r,a)) exp(-/^ ^ [[ ^ ^^Cr,©) - ik ^(r,a)) j exp(-/A:z)j 



[ni-6] 



where, for example, 
Therefore, 



a X 



} + ^{ 



V'A„(r,a))=j^ V'»I'(r,to) + ^»I'(r,to)-2jl |- ^(r,a))-A:' »I'(r,to)j exp(-/A:z) 



[111-7] 



[ni-8] 



and the parsed Helmholz equation ( without approximation) becomes 

Vf ^( r, co) - 2 jyt — ^ (r,a) ) + — ^ (r,a) ) = 

az az 



[ni-9] 



The paraxial approximation is precisely defined by the condition 



d_ 

dz 



lik — ^(r,a)) » — ^ ^(r,a)) 

oz 



[ ni-10 ] 
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which means that the longitudinal variation in the modulation function, ^(r,a) ) , changes 
very little in the wavelength associated with beam ~ i.e. In Ik . In this approximation, we 
neglect the third term and obtain the equation 



m{r,(ii)-lik — ^(r,(o) = 
oz 



[HI-ll] 



which is called the paraxial approximation of the wave equation. 

Solutions of the paraxial wave equation 

The Gaussian beam 

To inform or motivate our next step, we consider the paraxial approximation of a 
known solution of the Helmholz equation ~ i.e. a spherical wave 



exp 



{-ikr) exp(- ik ) 



exp 



■jx^ + y^ +z^ 




exp(-ifcz )exp 



-ik i^x + y^) 



2z 



[ III-12 ] 



Reflecting on the "quadratic" form of this approximate expression, it is reasonable to 
look for an axially symmetric solution of the paraxial wave equation in the following 
form ~ i.e. a Gaussian beam: 



^g(P, z, (o) = exp[-j P(z)] exp 



ik 

2^(7) 



[ III-13 ] 



Obvious the paraxial equation has the same mathematical form as the Schrodinger equation and, thus, all that is 
know about solutions of that equation may be directly appUed to understand issues in light propagation (or visa 
versa). 
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where =x^ . 

We test our conjecture by substituting the Gaussian beam function — i.e. Equation 
[ 111-13 ] ~ into the paraxial wave equation ~ i.e. Equation [ UI-l 1 ] ~ to wit 

expH Pfe)] v; jexp[-i^]- 2,-* A jexp[-,-p(,)] exp[-i^} =0 . [ 1.1-14 ] 



Executing the indicated operations, we obtain 



r ikp^il 2ik ^-Z .5 ikp^ 5 I A 



■lik 



-^"^^(^)+ .rTAf S^g(^) ^ = [IIl-15a] 



or simplifying 



^, 8 , lik ep^ [, a , J ^ 



[ III-15b ] 



Hence, for an arbitrary p this equation is separable into two parts ~ viz. 



— q{z)=l [in-15c] 

a z 



and 



^ , a , . 2ik 
2k—P{z) + --=0 
az q{z) 



a z ^(z) 



[lll-15d] 



which are satisfied by the simple solutions 

q{z) = Z+qo 



[ in-16a ] 
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and 

— P{z) = ^^ = -i — ln[z + qo] ^ P{z) = -iln[z+qo]. [ ni-16b ] 

oz Z+qQ oz 

On comparison with the paraxial approximation of a spherical wave — i.e Equation 
[ III- 12 ] ~ we may write q{z) in terms of a radius of curvature R (z) and a width 
w{z) -- viz. 

1 1 1 -i2 

+ TTT • [III-17] 



q{z) Z+qo R{z) k (z) ' 

To standardize the constants of integration we assume a plane wavefront at an 
arbitrary reference point z= ~ i.e. we take ^(0)=oo. Thus, 

^ =0 [III-18a] 



^(0) 

and => q^= = =i [ ni-18b ] 

k w [0) qQ 2 A 

where Lp = kw\0)/2 = 7U w (0)/A, is the critical Gaussian beam scaling parameter 
which is called variously the Fresnel length, the diffraction length, or the 
confocal parameter. In terms of this parameter, Equation [ 111-17 ] may be written 

— - — + r4^=^ = ^ [ni-19] 



q{z) R{z) kw\z) z + iL, z +L, 
Equating real and imaginary parts, we obtain 

^"774 ^ kw\z) ~ z+lI 
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R{z) = z[l^Llli] 
w^(z) = >v^(0)[l+z74] 

where Lp = 'Kw^ (0)/X 



[ III-20 ] 



Now since Equation [ III- 16b ] may be written 

P{z) = - i ln[z +qo\=-iln[z + iL,]=-i {in [z' + + i tan"' [^/z]}/ [ III- 16b' ] 



we may write 



exp[-jP(z)] = 



exp [-/ tan ' [Ip/z ]] exp [-/ tan ' [Lp / z]] 



to obtain the usual, officially approved form of the Gaussian Beam 
% ( p , z, CO ) = exp [ -i P (z )] expj^- j 



' 1 



. w{0) r „ -ht , ^^ ^ ikp^ 1 T p i 



[ ni-21 ] 
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The following kind of picture is sometimes found to be helpful in understanding the 
propagation of a Gaussian Beam (the bold curve depicts the spatial variation of the 
beam width and the light curve the beam curvature at particular points in space): 



4.T 



w(z) 
wCO) 



-3. 



-2. 



-1. 



1 . 



2. 




z it-VY'(O) 



Higher order Hermite-Gaussian beams 

In order to study the propagation of higher order beams, we substitute the following 
trial solution: 



%-g(P. z,a})=F{x, y, z) "^oip, z, co) 



[ III-22 ] 



into the paraxial wave equation - viz. Equation [ III- 11 ] — and obtain 
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F(x, y, z) Vf %(p, z, to) +2[v,F(x, y, z) ■ V, ^^(p, z, to)] +^g(p, z, to) F(:t; y, z) 



-2ik'¥^{p,z,(o)^ F{x,y,z)-2ik F{x,y,z)^ %(p,z,to)=0 
a z o z 



[ ni-23 ] 



Since the sum of the first and fifth terms already satisfies the paraxial wave equation. 
Equation [ 111-23 ] reduces to 



[ dw w~\ f g" [ dw w~\ g 



d^ 



— -\-2ik\— j^Y'' ^2ik\— \y 2kw — ;— =0 



\_dz q 



g 



\_dz q\ 



dz 



[ ni-24 ] 



From Equations [ 111-19 ] and [ 111-20 ] we see that 



dw w _w 
dz q~ R 



w —i2 
R kw 



i2 
k w 



so that the reduced equation ~ i.e. Equation [ III-24 ] ~ becomes 



g 



dz 



[ III-25 ] 



where ^ = xjw and c,=ylw . 



A Hermite polynomial of order n^^ has the following differential equation: 



^ H„ (t ) - 2t ^ H„ (t ) + 2n H„ (t ) = . 



dx 



dx 



[ ni-26 ] 



d" 

16 The Hermite polynomials have the generator H„(t ) = (-1)" exp(x ) exp(-X ^) . 

dx 
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With the simple change in variables T=^/2^ = ^/2x/>v and <3 = J[2c, =j2ylw 
Equation [ 111-25 ] may be written 



\\ dH ^ d^^ 1 \ d^g ^ dg^ 



- J 2 ~2X 1+ - I 2 I 

fL^^x dx \ g\_d(5 d(5 \ 



2a — I - 2 kw 



J4) 

dz 



= 



[ ni-27 ] 



Thus, it is apparent that we can write the functions ^^^^ ^(^~} Hermite 

polynomials ~ viz. 



f|-) = H„(T) = Hjy2^) and g|±) = H»=Hj^^ 



2 
w 



2d^ 



if we require that 2 A:>v = -2{n+m). Hence 

dz 

d^ _ (n+m) _ {n+m)Lp 
dz kw^ 2[Z^ + z^] 

or <i'{z) = -{n + m) atan{z/Lp). 

Finally we may write a general solution for the paraxial equation as 



o (P, z, (0) = A- ^ H„L/2 ^] hJ 72 1 



X exp[/ [n + m + 1] tan \z/ Lp)j exp 







r 1 


L 2i?(z)_ 


exp 


_ w^{z)_ 



[ III-28 ] 
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A Gallery of Hermite-Gaussian Field Distributions 
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Gaussian beam transformation matrices 

What we have shown above is that a given Hermite- Gaussian beam is essentially 
completely specified or defined by the complex function q{z). In propagating through an 
optical system, the beams are transformed by various optical components. The amazing 
fact is that the transformation produced by a given component follows a 
simple ABCD transformation law ~ viz. 



[ III-29 ] 




where A, B, C, D are the matrix elements found in our analysis of 
geometric optics!! 

To "prove" this, we argue by example. For example, the transformation through a 
uniform dielectric region is given by 

so that {A = 1; B = L; C = 0; D=1} and the transformation through a thin lens is given 
by 

J_- J__l 

?2 1l f 

SO that |a = 1; B=0; ^=-J'^ D=i|. 



Further "justification" of this transformation law may be found in terms of the so called 
"p/p' argument"— v/z 
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^r-i^in o f sphtrical 
wavt associo'ttd 
wi'th givtn inpu't ray 




■c^-fffmoi sphtrical 
wavt asf ocia'ttd wi'th 
givtn ou'tpu't ray 



From geometric optics and, in particular, Equation [ II-l 1 ], we may write 



Pout 
Pou, 



Pin 



C^+D 



or 



A^,=(AA^+B)(CAz,„+Dr 



which is identical to transformation equation [ 111-29 ] if we interpret q{z) as the wave 
optics generalization of A z = p/p ' . 
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IV. OPTICAL RESONATORS: 

Stability Criteria for Resonators and Periodic Optical 
structures by ray optic analysis 

Consider a prototypical periodic guiding lens system or an equivalent resonator. 

— d — — d 




m- 1 



m 



m+ 1 



Using the appropriate ABCD matrix with the indicated reference planes i^, we may write 

P... = Ap„+5p; [IV-la] 

and p:.i = Cp.+Dp; [IV-lb] 



1^ Where for reference , we see that 



'A B' 




"1 d~ 


1 0" 


"1 d~ 


1 0" 




C D 




1. 


-I//2 1. 


.0 1. 







l-d/f,-2 d/f, +d'lfj, 2d-d'/f, 
W2-V/1-I//2 1-^7/2 
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From the first equation we write 



Pm 



B 



and p;^i = 



B 



and substitute into Equation [ IV- lb ] to obtain 

P.. 2 - + ^] P..1 + [AD - BC] p„,, = 
The determinant of the coefficients \A D- BC\ = l so that 

P™.2-[^ + ^]P™+i + P™.i=0 • 



[ IV-2a ] 



[ IV-2a ] 



We see that 



A+D ^ d d d^ 
=|3 = 1 + — =-1 + 2 



/2 /i 2/,/, 



1-^ 
2X 



1- 



V 



2/2 



[IV-3] 



2 y 



Thus, stable ray propagation may characterized by bound solutions of the form 
Pm ~ Po exp(j ) which are possible if and only if 

exp(/(^) + exp(-/{^) = 2 cos{^ = A+Z)= 2|3 or cosc^ = |3 [ IV-4 ] 

Therefore propagation is stable ~ i.e. the rays are confined ~ when | P | < 1 so that 



0< 



1-- 



2/1 



1-- 



d 



2/, 



<1 



[IV-5] 



2 7 



Ray stability of rays in a periodic system may be usefully characterized in terms of the 

d d 

variables m, = 1 and = 1 as follows: 

2A 2/, 
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THE PARAXIAL WAVE EQUATION; 

Consider a Heraiite-Gaussian mode confined in an asymmetrical spherical cavity: 
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In order to sustain a resonant mode in such a cavity, the beam's radius of curvature must 
match each mirror's radius of curvature at the mirror's surface and, thus, the following 
conditions must hold (see Equation [ III-20 ] ): 



R,=z, + — and R2=-Z2-— [IV-5a] 



where d = Zi + Z2- Hence, we see that 



and Z2 = " . [ IV-5b ] 
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with a lot of algebra we can show that 



4 = 



[mj Ry -Mj^i] 



where now = l-d/R^ and Mj = ^-^dfR^ . 

For a symmetric resonator R^ = -R^ and u^= U2 



4 = 



_ d[2R-d] 
4 



and 



w(z,)=w(-Z2) = 



'dir 




-In J 


_d{R-d/2)_ 



IH 



[IV-6] 



[ IV-7a ] 



[ IV-7b ] 



For an asymmetric resonator, it can be shown with a bit more algebra that 



w(-Z2) 



'dX~\ 


if2 


-1; 


. n J 




Mj(l -Mj U2) 


'dX 










U2{l-U^U2) 



V4 



[ IV-8a ] 



[ IV-8b ] 



As a measure of the effect of resonator length and mirror radius on diffraction loss 
consider the ratio: 



=iRr^)\ 



[IV-9] 
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where >v_„f is beam width at the mirror for the confocal configuration ~ 

/. e. , when both mirrors have their focal points at the mid-point of the cavity 



2. 



1 . 





W f 7 ) 

1 w^ff^l^ 
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oonfigu 


ooa1 1 
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pUne-parallel 
configuration 

/ 

1 1 1 1 1 


concentric ■ 
configuration ■ 

\i 

1 1 1 1 1 ^ 



0.4 O.S 1.2 1.6 2. 

d 



Resonance Frequencies of the Optical Resonator; 

From Equation [ 111-28 ], we see the the "round-trip" cold resonance condition for a 
Hermite-Gaussian mode is given by 

kd-{n+m + l)[tan\zjLp ) + tan "'(z j/^f)] = [ IV- 10a ] 

where A'^ is an integer. In terms of frequency, the resonce condition is 



R. Victor Jones, February 10, 2000 



On Classical Electromagnetic Fields 



CO = -|A'^7U + (n+m +l)[tan '(zjlp)+tan '(z2/^f 



After much algebra, it can be shown that: 



to = ■^[^^ + (n + m + l)cos Mj] 
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[ IV-lOb ] 



[ IV-lOb ] 
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PLANE WAVE PROPAGATION IN A LINEAR, HOMOGENEOUS, ANISOTROPIC 
DIELECTRIC MEDIA ("CRYSTAL OPTICS"): 

Our objective here is to formulate a general approach to the subject of wave propagation in 
anisotropic dielectrics which makes use of ideas famihar from other branches of 
mathematical physics — viz. the "eigenvalue problem.". For reasons that will soon 
become abundantly clear, treatments of "crystal optics" focus on the behavior of the 
dielectric displacement vector, D(r, CO ) rather than on the electric field vector. For non- 
magnetic dielectrics the components of the dielectric displacement are usefully represented 
as the Cartesian coordinates coordinates of figure called the "ellipsoid of wave normals," 
the "optical indicatix," the "index ellipsoid" or the "reciprocal ellipsoid." 



Since the stored electrical energy is given by 

f/, = (l/2)DE=(l/2)D-8"'-D 

we can write in the principle axis system 



D. 



V 



D. 



V 



+ ■ 



D_ 



V 



= 1 



where {e^ , £,,,£,} are principle axis values of the dielectric 
constant tensor. 




Kaiser S. Kunz in 1977 presented a similar treatment of this problem in a paper enti tiled "Treatment of optical 
propagation in crystals using projection dyadics," Am. J. Phys., Vol. 45, 1977, pp. 267-269. 

Perhaps the most authoritative treatment of "Crystal Optics" is found in Max Born and Emil Wolf, Principle of 

Optics, Pergamon Press (1986), Chapter 14. r. victor jones, December 19, 2000 
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Now let us first combine Equations [ I-8a ] and [ I-8b ] to obtain a generalized Helmholz 
equation for a homogeneous anisotropic dielectric 

Vx(vx[r'((o)-D(r,(o)]) = (oV« iXr,co) [V-1 ] 

In general, we would hope to be able to find a set of eigenmodes of the homogeneous 
problem that would satisfy the scalar eigenequation 

V'D"'\r,03) + nl^^kl D'"\r,a))=0 [V-2] 

where n^^^ is the effective index of refraction of the a-th eigenmode. We can be quite 

specific for the case of plane wave eigenmodes where we suppose that all fields have a 
spatial dependence exp(+ / k r) . From Equation [ I-8c ] we see that D(r, co) must, in 
general, be orthogonal to k so that we may write 

D (k, CO) = D^''(k, CO) t'-'Xk, co) + D'''(k, co) t^'^(k, co) [ V-3 ] 

where the t*°' (k)'s are polarization unit vectors which are orthogonal to k , the unit vector 
parallel to k . Since the components of 8 in the general case may be complex ~ e.g., in the 
case of magneto-optical media ~ the eigenmodes may be polarized along "complex 
directions" ~ e.g., "screw axes" defining the sense of circular polarization ~ and we must 
use great care in all vector manipulations. Thus, we define a set of adjoint or conjugate unit 
vectors by means of the relationships 

k-k=l; k-t('"(k)=0; t(<"(k)-k = 0; ^''^(k)- t<''''(k)= 5„„, [ V-4 ] 

Thus, with k = k k = k* k (so that k k = |k|^ ) the representation for D(k) in Equation 
[ V-3 ] automatically satisfies Equation. [ I-8c ]. A key problem is the representation of the 
ubiquitous vector operation 
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3 (r)-F(r)= V x(v xF(r)) = W •F(r)-V'F(r) [ V-5 ] 

which appears in Equation [ V-1 ]. For plane wave representation of the operator becomes 

3(k) = (k-k)i-kk = |k|'{i-kk}. [V-6] 

Since 

1 = k k + t''\k) V\k) + t*" (k) t*" (k) [ V-7 ] 

the plane wave representation of the operator simplifies to 

3 (k) = |k|' {t^'\k) t^'\k) + t'"' (k) t'"' (k)} [ V-8] 

Using this representation of the 3 (k) operator and the representation for D(k) in Equation. 
[ V-1 ], the generalized Heknholz equation ~ i.e. Equation [ V-1 ] ~ can be written 

ikf {D<"(k) v\k) ■ r • t^'^(k)] + (k) (k) [t (k) ■ ■ t^'^(k)] + 

D (k) t^''( k) [ t *'\k) • 8"' • t (k)] + D^'' (k) t^'' (k) [t*'' (k) • e"' • t''' (k)]} [ V-9 ] 
= 0) Vo{D''\k) P'(k) + D*^> (k) F(k) } 

Crucial point ; To obtain an eigenvalue equation we need to choose the eigenvectors 
t^°^(k) so that 

t^'\k) • e~' • t^" (k) = t^'' (k) • r' • t*"(k) = [ V-10 ] 

If we can find eigenvectors defined in this way, then Equation [ 1-9 ] becomes an 
eigenvalue equation with eigenvalues - i.e., the inverse refractive indices - given by 

[n^(k)]"' = k^/|k|'=8o P\k)-r'-F\k) [ V-lla] 
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[n,(k)]" =kl/\kf = e,V'\k)-r''e\k) [V-llb] 

These results ~ i.e., Equations [ I-IO ] and [ I-ll ] ~ are a complete formal solution of the 
problem. However, they are difficult to apply in the general case and an addition 
relationship ~ viz., the Fresnel equation of wave normals ~ is found to be extremely useful 
as the starting point for actual computations. For plane waves. Equation [ V-1 ] can be 
rewritten as 

[|k|' i - (koVeo) e] -£(1^) = k k •E(k) [ V-12 ] 

Multiplying this equation through by k • pk|^ i - (kg/eo) 8 j • we obtain 

k £(1^)= {k- [ikf i -(kl/E,)eY ■ k} k- E(k) [ V-13a ] 

or l=k-[|kf i-(ko/eo)8]"'-k [V-13b] 

From this result we may develop two important relationships. Using the principal axes 
coordinates of the dielectric tensor, we can write 



[ikp 1 - (koVeo) e] = - (ko/eo) \ a, [ V-14a ] 

X= a 

and [|k|' 1 - (klle,) e]"' = f^(\kf - (k^ej E(,))"' a, [ V-14b ] 



Therefore, Equation [ I-14b ] becomes 



(k.aj(k.aj (k.a,)(k.aj (k.a)(k.aj 1 

22 22 22 2 >- J 



n —n^ n —n^ n —n^ 
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where nl = £(x)/^o ■ Since a^^ = 1 (or ^(k • )( '^)-^ ) we may also write 

X=a X=a 

Equation [ V-15a ] as 

(k-aj{k-aj {k-a,){k-a,) {k-aj{k.aj 



1 1 



1 1 



+ ■ 



1 1 



= 



[ V-15b] 



2 2 

n n 



2 2 

n Hu 



2 2 

n n 



This latter expression is the famous Fresnel equation of wave normals. 

APPLICATIONS OF THE FORMAL SOLUTION: 
Uniaxial Dielectric Crystals : 

For an optical material with uniaxial symmetry, the inverse dielectric tensor in the 

principal axes system must have the form 



[V-16] 



Thus, Equation [ I-19b ] becomes 



j__J_ 

2 ~ 2 



n n 



1_J_ 

2~ 2 

n «|| 



sin 1^ + 



j__J_ 

2 ~ 2 



n n 



cos 13 r = o 



[V-17] 



so that 



=nl\ and =n'l cos^i^ + nn"^ sin^-d 



[V-18] 



[^l H n 

2" In its commonly used form, the Fresnel equation becomes ~ + ~ ~[ ~ ~ ^ 

2 2 2 2 2 2 

n n n 

21 In this instance there is no need to trouble ourselves about conjugate unit vectors. 

R. Victor Jones, December 19, 2000 



On Classical Electromagnetic Fields 



Page 47 



where k- 3^ = 00813^. The subscript "o" identifies the "ordinary" mode and the 
subscript "e" the "extraordinary" mode. These results are usually plotted as follows: 




where the intersections of the k vector yield the "ordinary" and "extraordinary" 
velocities of propagation for a given k . 

Further, if we take 

k = sin(p sini3 + cos(p sini3 a ^ + cosi3 a^ [ V-19 
it is a bagatelle to show that 

t*°'(k)= r^'Ck) = cos(p a^ - sincp a, [ V-20a ] 

t*'' (k) = t (k) = cos-d (coscp a^ + sincp a J -siniJ a, [ V-20b ] 

and that these quations are consistent with Equations [ V-10 ] and [ V-11 ]. 
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See, for example. Section 82 in L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, 
Pergamon Press (1960). 
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i+tY(a„a,-a,aj]8-^ [V-21] 



If we introduce the conjugate principal axes 

a.=a_=(l/^)(a„+ta,) 

5_ = a, =(l/^)(a,-taj [V-22] 

we obtain the dielectric dyadic in the so called normal form -- viz. 

= [(1 - Y ) a, a, + (1 +Y ) S_ a_ + 5„ S„ ] e"^ [ V-23 ] 

Again from Equation [ V-15b ] it is trivial to show that 

=[1±Y cosT^]8"' [V-24] 

Using the resolution of k as given in Equation [ V-4 ] we may show that 

P (k) = (^^ ^ exp(- 1 (p) [ 1 + cos 1^] a^ - exp(i (p) [1 ± cos a_ ± ^2 sini^ a,, } [ V-25a ] 

and 

t-(k)= [^^j{-exp(i(p)[l+cosi^] a+ +exp(-i(p)[l±cosi^] a_ + >72 sini^ a,,} [ V-25b ] 

ENERGY FLOW IN ANISOTROPIC MEDIA 

As previously noted, the content of Equations [ I-IO ] and [ I-ll ] represents in some 
sense a complete formal solution of the wave propagation problem. However, from a 
practical point of view it is essential to consider how energy propagates in anisotropic 
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media. To that end, we note that the time averaged Poynting vector associated with a 
given plane wave- viz. 



S (k, (o) = ^ E (k, (o) X H* (k, (o) 



[V-26] 



propagates in a direction s (k, (oj — conventionally designated the "ray" direction 
which is orthogonal to both E (k, co) and 11 (k, O)) as shown below. 




The time averaged total stored energy is given by 
1 r - /r 



t/(k, co) = - [e (k, (o) ■ D* (k, co) + B (k, co) ■ H* (k, co)] 



= ^k E(k,co)xH* (k,co) = ^'k S(k,(o) 



[ V-27 ] 
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and, thus, we see that the "ray" or "energy flow" velocity, v^^y , for a given s (k, co) i 
given by 



IS 



1 I ^ \ 1 ^ \ 
— = — k s k, CO = k s k, CO 



[ V-28 ] 

' ray ™ ' phase 

We write the time averaged Poynting vector associated with a given eigenmode as 



S"" (k) = i E"" (k) X H'-i* (k) = E"" (k) X [k X E"'* (k)' 



[ V-29 ] 



where E*"' (k) is the electric field associated with the eigenmode. Using Equations 
[ V-3 ], [ V-10 ] and [ V-U ] this field can be expressed as 

E<"' (k) = {[t'"\k) • 8"' • t<"\k)] V'ik) + [k-r' ■ t'"\k)] k} D'"'(k) 



r ^ ^.2 {t^'^kk) + Eo [«(a)(k)]'[k- e"' • t'^^k)] k} D'^^k) 
[«(c)(k)J 



[ V-30 ] 



^ {t'">(k) + g^">(k)k}D'°\k) 



where g<"\k)= eJ.^.^Ck)] [k.rM<^>(k)] = ^1^^^^^^ 
parameterization, the modal Pojmting vector can be expressed as 



Using this 



[V-31] 
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and the associated ray vector as 

k-t<'"(k)^^'"(k) k-P'(k)^^'''(k) 



Jk-t<''>(k)g^<'>(k)| Ji+ig^^'ck)!' 



If we take k and t'^^k) as reference directions, the ray vector, s*"' (k), lies in the 
plane containing k and t'^'''(k) at an angle T|*'"(k)= -tan'^fg^^^k)] with respect to the 
direction k. 

For an optical material with uniaxial symmetry, we may use Equations [ V-16 ], 
[ V-19 ] and [ V-20 ] to evaluate g'''^(k). In particular, we may easily see that for the 
ordinary mode 

tami*°\k) = -g^°'(k) = [ V-33a ] 

or s^°' (k) = k and for the extraordinary mode 



sin-dcosi3^ (£± 

e;^' cos^-d + e,;^ sin^iJ' 



tanri^^^k) =-/\k) = ^y.]:''^,^^^ . [ V-33b ] 
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VI. DESCRIPTIONS OF POLARIZED LIGHT22 

Consider a totally coherent wave propagating in the positive direction 



E, = E,' cos(a)? +6 J = E° [exp[/ (to? + 5 J] + c.c] 
E^ = E^° cos (co r + 5 J = e; [exp [j ( cor + 5 J ] + c.c] 



[ Vl-la ] 
[VI- lb] 



For later reference, we note that the full (non-normalized) Jones vector representation ^3 
of this field is given by 



E° exp id 



[VI-2] 



We can be easily shown that 



cosa)r=[sin(5^-5j] j^|s-sin5^ 



smco? = 



rE 



= [sin(5^-5j] l^^cosS 



E 1 

■^sin5,j 



E 1 
■^cos6,j 



[ VI-3a ] 



[ VI-3b ] 



The best references on this subject are the following: I.) William A. ShurcUff, Polarized Light: Production 

and Use, Harvard University Press (1962); 2.) D. Clarke and J.F. Grainger, Polarized Light and Optical 
Measurement, Pergamon Press (1971); 3.) Max Born and Emil Wolf, Principles of Optics, Pergamon Press 
(Particularly Section 1.4). 

R. Clarke Jones, "New calculus for the treatment of optical systems. I. Description and discussion of calculus," 
/ Opt. Soc. Amer. 31, 488 (1941). To quote Shurcliff, 

"The Jones vector, . . . describes a polarized beam with the maximum algebraic brevity, and is 
eminently suited to the solving of problems involving beams whose phase relations phase 
relations are important. ..." 
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Thus the locus of time sequence of fields in a plane perpendicular to the direction of 
propagtion follows an ellipse - viz 




The 'plane of polarization' vith propagation out of the plane 

I 

For an electric field vector "seen" to be rotating in a clockwise direction by an observer 
receiving the radiation (/.e., n> (5^, -5j> 0), the polarization is said to he right-handed. 

For rotation in the anticlockwise sense (/. e., 0> (5^, - 5 J >- 71 ), the polarization is said to 
be right-handed. Although this is a complete decription of the coherent field, it is not a 
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convenient one. It is useful to transform this equation to its principal axes form by the 
following transformation 



= cos(t) -E^ sin^ 



Ey = E^ cos(|) + E^ sin^ 



[ VI-5a ] 
[ VI-5b ] 



where we choose (|) so that 



^E. 



+ 



^E^ 



= 1 



E^ = a cos(o)?+6o) 
E^ =±Z? sin (cor +5 o) 



[ VI-6a ] 

[ VI-6b ] 
[ VI-6c ] 



Where the upper and lower signs are for, respectively, right-handed and left-handed 
polarizations. After quite a bit of algebraic manipulation, we find a more elegant and 
convenient description of polarization in terms of the following set of relationships 



tan2(t) = 



aZ, = (E:)(E;)sin(5,-6j=(E:)(E;)sin6 
2(E«)(E°)cos(5^-5j 2(e1')(e;')cos5 



2ab 2(E;)(E^»)sin(5^-5j _ 2(E,»)(E';)sinS 



2 , 1,2 

a + b 



[ VI-7a ] 
[ VI-7b ] 

[ VI-7c ] 
[ VI-7d ] 
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These relations are simplified if we introduce the two auxiliary angles 

tana=EyE°, [VI-8a] 

and tane =± bja, [ VI-8b ] 

b - a 

(Note: — — specifies the, so called, ellipticity of the vibrational ellipse.). In terms of 
these auxiliary, we may then write 

[K)'AK)'=a^b' [VI-9a] 

2tana „ / \ „ 

tan 26 = :; — cos6 =( tan 2a I cos5 [ VI-9b ] 

^ 1-tan a ^ ' 

sin2B =( sin 2a ) sin 5 [ VI-9c ] 

Probably, the most powerful representation of polarization is found in the famous Stokes vector or 
parameters ~ v/z.^^ 

/=(E°f + (e;)' [Vl-lOa] 

M=(E°)'-(E';f [Vl-lOb] 

C = 2 (E° ) (e;) cos (6, - 6 J = 2 (e") (eJ) cos 5 [ VI-lOc ] 

S = 2 [El ) (e;) sin (6, - 6 J = 2 (E^ (eJ ) sin 5 [ Vl-lOd ] 
so that a polarized field may be represented by the vector 



G. G. Stokes, "On the composition and resolution of streams of polarized light from different sources," Trans. 
Cambridge Phil. Soc. 9, 399 (1852). 
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/ 

M 
C 
S 



or its transpose 



{l M C S} 



[ Vl-lla] 



[ Vl-llb] 



These vector components give a complete geometric description of the vibrational ellipse ~ viz. 

I Size 



C/M = tan2<|) 

N//=sin2e = / ,2 
l + {b/a) 

Sign of 5 



Azimuth 
Shape 

Handedness 



As can see from Equations [ VI- 10 ], for a completely polarized field only three of these 
parameters or vector components are independent, since 



[VI-9] 



and {M, C, S} can be interpreted as the Cartesian coordinates of a sphere of radius I ~ the 
Poincare sphere. ^5 Thus, from Equations [ VI-8 ] and [ VI- 10 ] we obtain the coordinates 
of the Poincare sphere or representation as 



H. Poincare, Theorie Mathematique de la Lumiere, Vol. 2, (1892) Chap. 12. 
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/'=M' + C' + 5' [Vl-lOa] 

M=/cos2e cos2(t) [Vl-lOb] 

C=/cos2esin2(^ [Vl-lOb] 

5=/sin2e [VI- 10b] 



R. Victor Jones, February 17, 2000 



On Classical Electromagnetic Fields 



Page 56 



LHCP 




RHCP 
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[VIM] 



VII. NONLINEAR Optics -- Classical Picture: 

An Extended Phenomenolqgical model of polarization; 

As an introduction to the subject of nonlinear optical phenomena, we write, in the spirit of 
Equation [ 1-4 ], the most general form of higher order terms in the phenomenological 
electric field expansion of the polarization density (which may then be inserted in Equations 
[1-3 ]) as 

K''\^^t)= EoLJ I I \d%dt,di,dt,x'^l^{i-i, t- t{,i-%, ?-?2)Ep(?,^)E^(f2, t^] 

+^0 L I I I I I I d% dt^ dr^ dt^ xl^^p^g (f- r,, t- t{, f- t- t^;r- fj, t- 

xEp(fi,/i)E^(f2, t^] Egffj. r3) + ---. 

The wave vector and frequency dependent second and third order susceptibilities 

are then defined as 

X^j^^kj, (Oi;k2, 0)2]= J J J J^i?Rii/Xii/R2<iX2exp[-j kj-Rj exp[+j(Oi xj [VII-2a] 

Ri R2 '^i ^2 

X exp [-i k 2 • R 2] exp [+i (O2 X2] X Sy [% '^i; ^ 2- ] 

and 

Xapys k2,C02; 63,0)3]= J J J J J J dti^dx^ dR2 d'^2 di^ exp[-j kj-Rj exp[+i(0^ xj 

Ri R2 R3 ^1 ^2 1^ VII-2b ] 

X exp [-/ kj • R2] exp[+ia)2 X2] exp 63 • R3] exp [+/CO3 X3] 

xxSy8(Ri.'Ci;R2. X2;R3.X3 
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Thus, we may write quite generally 



26 



Pj^""(r, = Efl 12 1 1 J J JcOj exp| / ( kj+ k^] ■ rj exp[-/ (co^ +(0^) ■ t] 

X% apy [K ^2 ^2) Ep (kpCO,] E^(k,, (0,) 

+eo^ 1 1 J J J Ji/kji/cOj dk^ do)^ dk^do)^ exp|/ (kj + k^ + kjj fj exp[-/ (cOj+cOj+cOj) 

PtS ^ kj ^ (0, (02 

x5Capy8 (ki.c0i;k2. (>i2,K cOjjEp (ki, co) E,(k2.co2)E5(k3,co3) + ---. 

A Simple classical model of nonlinear optical response 

A simple Lorentz-Dude model is often used in the literature as a valuable guide to the 
understanding of the frequency behavior of the nonlinear dielectric response.^^ We assume 
that the potential energy of a one-dimensional nonlinear (anharmonic) oscillator may be 
written 

V{x) =V^{x) + V,{x) + V,{x) + --- = ^M(olx^ + ^Ma x^ + ^Mb x' + --- [ VII-4 ] 

(See figures on next page) 

Thus, the equation of motion of a particle moving in that potential becomes 

x+Fx+dilx + ax^ + b x^ + ---= — E(t) [ VII-5 ] 



[ VII-3 ] 



% pjP^ must vanish for any material that is invariant under inversion, since both P and E are vectors, and are 

(3) 

thus odd under inversion symmetry. Note also that % apys ^ given material has the same transformation 

(2) (3) 

properties the elastic constants of that material. The nonzero elements of Xaf,y XapyS ^'^^ various crystal 
symmetries are compiled in Y.R. Shen, The Principles of Nonlinear Optics (Wiley, New York, 1984). 

See, for example, N. Bloembergen, Nonlinear Optics (The Advance Book Program), Addison- Wesley (1992), 
ISBN 0-201-57868-9. 
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We can analyze the response of the oscillator by expanding the displacement in powers of 
the electric field E[t) — viz. 

it) = x''\ t) + /'\t) + x''\t) + ■■■ [ VII-6 ] 

where x^"\t) is proportional to the nth power of the field E (?) . Inserting this expression 
into Equation [ VII-5 ] and equating like powers of E(?), we obtain the following 
hierarchy of equations: 

x%)+rx^'\t)+&l/'Xt) = {q/M) E{t) [ VII-7a] 

x''\t)+rx''\t)+(iil /'\t) + a[x'\t)f = [ VII-7b ] 

x''\t)+rx''\t)+(olx''\t) + 2a x'\t)/'\t) + b[x''Xt)]' = [ VII-7c ] 
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In the frequency domain we see that 

-co' /'\(£)) - m Tjc^'^co) +03 ^ jc'"((o) = {q/M) E (co) [ VII-8a ] 

or xjoi) = (^/M) E(a)) [(o^-co^- /cor] = (^/M) E((o) ^((Oo,(o;r) [VII-8b] 

where ir(M,v; w)= [m^- v' - / v w] .^^ Thus, knowing x^^'(a)), we may then consider 
x^'^(r) a driving term in Equation [ VII-7b ] ~ viz. 

x''\t)+rx''\t)+(olx'y) =-a [x''\t)f [ VII-9 ] 

Therefore, in the frequency domain 

x''\oi) = -a !F{(0,M I) I Jco' x''\(o) x''\(0- (o) ^ ^"'^^^ ^ 

which, in view of Equation [ 1 l-8b ], becomes 
x'^^\(i)) = -a{q/Mf ir[cOo,CO;r] 



X I Jco' [cOq, co'; r] [cOq , (co- co') ; r] E (co') E (co- co') 



.[ Vn-lOb ] 



We may treat the third order terms in a similar manner. We write Equation [ Vn-7c ] as 

x''\t)+rx''\t)+(Dlx''\t)=-2a x'\t)x''\t)- b[xy)\' [VII-U] 



Note that near resonance 



where (D is the so called complex Lorentzian. 
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In the frequency domain 

x^'\(ii) = -2a ir(a)o,a); I) J dd)' x^^\(xi- to") . VII-12 ] 

- & (tOg ,to; r) J J dca' da" jc'^^co') jc'^^co") jc^^^co- co'- co") 

Using Equations [ VII-8b ] and [ VII- 10b ], we obtain 

x^'\(o) = 2a' [q/Mf !F[oio,(o;r] J J d(o' dco" !F[(Oo,{(0- a)');r] iFK,©'; T] 

X ^[a)o,(o"; rl ^[(Ooio)'- (o"); rl E (co- co') E(co") E(co'- co") 

ff [VII-13] 
-b{q/My JK,co;r] J J ^ico'^io)":r[o)o,co';r] rK,o)";I] 

X ir[a)o,(a)- CO - co"); r]E(co') E(co") E (co- co - co") 

In particular, for an input 

E (?) = |Ej| cos(o3i?+(p i) + |E2|cos(c02?+(p 2) 

= i [Ej exp(- + E* exp(+/a)i?)] + ^ [^2 exp(-/C02?) + E* exp(+/C02?)] 

a. Second harmonic generation is due to the terms 

x^''(2cOi) =-a{q/2Mf ^ [co^ , 2c0i; T] ^'[cOo,co,;r] E'(cOi) [ VII-15a ] 

and 

x^'^(2co2) =-a [q/iuf :r[cOo,2c02;r] ^'[cOo,c02;r] E2(co2) [ Vll-15b ] 

b. Sum frequency generation is due to the term 
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x*^^(to,+a),) = -a (q/2MY irfa)n,(a),+a),);rl 

V 1 2/ J I oy I 2j i ^ VII-16a ] 

X :r[a)o,a)i; T] J^[(>}q,()}^; T] E/coJ £2(0)2) 

c. Difference frequency generation is due to the terai 

x^^Vto.-co,) = -a iq/lMf ^[a)„,(co,-a),);rl 

X ir[a)o,a)i; T] *[a)o'f02; T] Ei(a)i) E*(co2) 

d. Optical rectification or dc generation is due to the tenns 

x''\0) =- a {q/2Mf ^[(>i„0; T] \9^[oi„Oi,;r]( |E,((ojf [ Vll-lVa ] 

and 

;c*''(0) =-a{q/2M)' ^K,0;r] {j^ia^^,; r]f \E,{(>i,)( [ VII-17b ] 

e. Third harmonic generation is due to the tenns 



X [cOo , 3 coi; T] r [tOo , <Si{, T] Ef (cOi) 

and 

x^^\3(ii^ = {qlMf {2fl2 ^[C0o,2co2;r]-Z?} 

xirK,3a)2;r] iPK,a)2;r] £^(0)2) 

f . Intensity dependent propagation is due to the tenns 

x(3>((oJ = [qjlMf [la^ ^[a)o,2a)i;r]- 

x\j [0)0,0)1; r]f r [o)o,0)i; r] |Ei(o)i)fEi(o)i) 



[ VII-18a] 



[ VII- 18b ] 



[ VII-19a ] 
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and 

x<3)(a)2) = {q/2My {2a' ^F^.l^^ir]- 

2 [VII- 19b] 

x|i^K'0)2;r]|Y'[coo,co2;r] |Ei(co2)| eIo,^) 

g. Raman generation (inelastic scattering) involves tenns like 

20)1-0)2)= 2(^/M)^{2a' :r[o)(,,(o)i-0)2);r] + a' 9^[(i3Q,2(i3{,T\- b} 

xiF[o)o,(2o)i-0)2);r]i?^[o)o,0)i; T] f [0)0,0)2; T] Ef(o)i) E*^) 

Note that, according to this simple anhannonic oscillator model, Raman generation may be 
enhanced by a resonance at a frequency {(O^ - (i>2) ! Also note that, for this model (see 

Equation [ Vll-lSa ] above), the ratio 

is a constant independent of frequency ! This observation is consistent with the famous 
empirical Miller Rule which declares that the ratio 



Xg.((03=tOl+(02) 

x1!\q)3) x^N 



A = — ' ■ — r VII-22 1 



has only a weak dispersion and is almost a constant for a wide range of materials! 

Second harmonic Generation -- Perturbation Analysis 
We may write the nonlinear , macroscopic Maxwell equations in the form 

V'E(r,0))- V[V •E(r,0))] + (o)Yc'8o)l E(r,0))=-^ioO)' P*'"''(r,0)) [ VII-23 ] 



[ VII-20 ] 
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V {i ■ E(f ,0))] = -V •P^'^^f ,0)) [ VII-24 ] 

Suppose that we have an input driving or pump field 

E(r,C0i)= e*" E(cOi) Qxp(i{k^z-(iiit)) + c.c. [ VII-25 ] 

Then the components of the polarization with frequency 2c0j are given by 

P^'^\z,(02= 2a)i)= E'(a)i) |'":e*" e^'' exp[/(2;^z- 2K0,t)] + c.c. [ VH- 26 ] 

It is a convenience to resolve this p^'^''^ and the resultant second harmonic field into 
longitudinal and transverse components ~ viz. 

P*^"'(z,(0,)= zP;^^'(z,03,) + Pf"'(z,(0,) . [ VII-27a] 

and 

E(z,a)2)= zE,|(z,a)2) + E^(z,a)2) [ VII-27b] 

Therefore, we may write the cOj or second harmonic components of Equations [ VII-23 ] 
and [ Vn-24 ] 

|1 E,(z,Q),) + ^ ^ E,(z,G),) to^ P,'^^'(z,a),) 
az c Eq 

° 1 [Vn-28a] 

+R ^ E„(z,a),)+Hoa)^Pr(z,a),) z= 

I ^ ^0 J 

and — |e|,((02) E|,(z, 0)2) + P|r'(z, 0)2)1=0 [VII-28b] 

oz 

To satisfy these equations two conditions must hold ~ viz. 
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E,(z,a)2)=--7|-Tlir(z,03,) 



[ VII-29a ] 



and 



[ VII-29b ] 



where 



^2= — x/e±K)/eo= 2cOi J|ioe^(2(0i) 



[ VII-30 ] 



We now write Ej^z,(i>2)- ^'^^ £(2,032) exp( i k2z) where £(2,0)2) is a slowly varying 
function ofz ~ viz. 



^E^(z,C02) = -e''' 1^2' £(2,0)2)- 2/^2 ^ 2:(2>«2)- ^ 2:(2>«2)|exp( ^^2^) 

£(2,0)2)- 2/ /:2^ £(z,0)2)|exp( 1^22) 



[ VII-31 ] 



so that Equation [ VII-29b ] may be written 



^£(2,0)2)= — lioCO^ e<2':P^<^^)(z,0)2)exp(-//:2z) 



dz 



2k. 



_ tO; Tj2/ X -(2). =(2) -(1) -(1) 



[ VII-32 ] 



2L 



E>i)e -I -e e exp(iMz) 



where Ak= 2^- ^2 = — ^ ['^fh i'^) ~ "sh (2o)i)] 
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If we assume that the driving field stays constant, we can directly integrate 
Equation [ VII-32 ] to obtain the exceedingly famous and important equation for 

the spatial variation of the second harmonic field ~ viz. 



£(2,0)2) = z 



'"M-O -(2). i(2).-(l) -{1) 

4L 



^[2), ^[2),^(i) ^(ij exp(/M z/2) ^E'(tOi) sinc[Mz/2] .[ VII-33 ] 



MAKER'' FRINGES 




-40 -20 20 

ANGULAR ROTATION IN DEGREES 



40 



LASER BEAM 



SECOND 
HARMONIC 




z = d / COS e 



Second harmonic Generation - Coupled Wave Analysis 

When the process of second harmonic generation takes place under conditions of perfect 
phase match — i.e. k2=2ki -- the perturbation result breaks down if the path is 
sufficiently long. Under these circumstances the pump beam will be depleted as the second 
harmonic grows and a solution of Equation [ VII-32 ] must take into account the spatial 
variation of E' (tOj) . To that end we assume perfect phase matching 



P.D. Maker, R.W. Terhune, N. Nisenoff, and CM. Savage, Phys. Rev. Lett., 8, 21 (1965). 
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n,(2co) 



n,(2co) 




Tl/2 



and rewrite Equation [ Vn-32 ] as 



OZ 



[ VII-34 ] 



Of course, as the second harmonic grows Equation [ VII-3 ] tells us that a nonlinear 
polarization at the pump frequency is generated — viz. 



P^^^'(z,a)0=E{z,2(o,)E^(z,a), 



[ VII-35 ] 



Repeating the arguments associated with Equations [ VII-28 ] through [ VII-3 1 ] we may 
write 



az 2ky 



[ VII-36 ] 
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Combining these two equations, we find an equation governing the spatial evolution of 
E(z,a)i) - viz. 



E(z,2a)i)E*(z,a)i) . [VII-37] 



Equations [ VII-34 ] and [ ¥11-37 ] are then the coupled differential equation which 
describe the coupling of the first and second harmonic fields. Handling all the "vectomess" 
in these two equation would obscure important issues. Thus, we consider a pair of 
somewhat less complicated equations which incorporate the essence of the problem ~ viz. 



OZ 



[ VII-38a ] 



|-E(z,coO=%^x'^'E(z,2cojE*(z,co.) 
OZ Aq 



[ VII-38b ] 



To solve these equation we first write E (z,©;) = E (oJi) fi(z) and E (2,0)2) = E ((Oj) f2(z) 
with the boundary conditions fi(0) = 1 and fjCO) = 0. Thus, the coupled equations reduce 
to the dimensionless form 



[ VII-39a ] 



and 



f^fi^)=imm 



[ VII-39b ] 



where ^= z/L^ and lJ= (of%*^' E((ojjy^. We next separate fi(^) and f2(^) into 
phase and amplitude parts as 
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l2(^)=«i,2(^)exp[/(^i,,(^)] [VII-40] 
and substitute into Equations [ VII-39a ] and [ Vn-39b ] ~ viz. 

-^«,(^) + /«,(^)^(t),(^ = ^>^(^)exp{/[2^,(^)-^,(^)]} [VII-41a] 



j-^u,{^) + iu,{^) = i a^(^)exp{-/[2(^,(^)-^,(^)]} [ VII-41b ] 

Equating real and imaginary parts of these equations, we find 

^ u,{^) = -u^i^) sin [2^,(^) ,(^)] [ VII-42a ] 

f^U^) = «i (^) cos [2^,(^) ,(^)] [ VII-42b ] 

^ u,{^) = u,{^) u,{^) sin [2^,(^) ,(^)] [ VII-42C ] 

^^^^) = cos[2(^j(^)-(^,(^)] [ VII-42d ] 

Since the phase enters only in the combination \|^^) = 2(t)i(^) -<t> il^) these four equations 
reduce to three ~ viz. 

j-^u,{^)=-u,^)sm^^f{^) [VII-43a] 

u,{^) = u,{^) u,{^) sin V|/(^) [ VII-43b ] 
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Combining these three equations, we obtain 



which is, obviously, equivalent to 

^ In [ u,{^) uf{^) cos \|/(^) ] = [ VII-44b ] 

or 

M2(^) cos\|/(£,) = const. [ VII-44C ] 

Since MjI^) =^ as the "const' must be zero and, hence, \|/(^) must be %/2 for all 
^. Thus, the original four coupled equations now reduce to two ~ viz. 

j-^u,{^)= u,\?,) [VII-45a] 

and 

j-^u,{?,)=-u,{^)u,{^) [VII-45b] 

Combining these equations, we obtain 

ui^) ^ u,{^) + u,{^) ^ u,{^) = ^ + 4^)] = [ VII-46a ] 



or 
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uf{^) + ul{^)=l [VII-46b] 

which is an assertion of the principle of energy conservation. Taking this equation 
together with Equation [ VII-45a ] we see that 

-^ = d^ [VII-47a] 



or M2(^)=tanh^ [ VII-47b ] 

and JlnM.(a=-tanh(a => m,(^)= — ^— [VII-47c] 

cosh(g 

Finally, returning to the original variables, we see that 

E{z,(xi^) = i E(a)i) tanh(z/Lj [ VII-48a ] 

/ N E(co,) 

and E{z,(0,)= J ][ . [VII-48b] 

cosh(z/LJ 
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VIII. Guided Waves in Planar Structures 

Characteristics of Plane Wave Solutions; 

For the record, let us restate the frequency domain, macroscopic Maxwell's equations which are 
valid in the high frequency or optical regime for a linear, local, isotropic medium — viz. 

VxE(f ,co) = V X8~'(f ,(o) D(r,co)=-/coB(r,co)=-/co|j.oH(r,co) [VIII- la] 

VxB(r,a)) =iIq VxH(r,a))= ix^ j(f,a)) + /a)|i.o8(f,a)) • E(f,a)) 

-w^ ^ ^ [VIII- lb] 
=|io J(r,a)j + ia)|ioD[r,a)j 

V e(r,a)) E(r,a)) = V D (r,a)) = p (r,co) [ VIH- Ic ] 

V B (r,©) =|i V •H(r,a)) = [ viH- Id ] 

Further, in regions free o f explicit sources of current and charge we may write 

V xE(f ,co) =- / co|ioH(r,(o) [ yill- 2a ] 

V xH(f,a)) = / a)8,ff(f,a)) • E{m) ^ ] 

V e(r,to) E (f ,0)) = [ vill- 2c ] 

V •H(r,a)) = [ viH- 2d ] 

where 8^(f(r,(o) =8 eff (? .w) - i o (f ,(o) /o) . In this set of lectures it is our intention to explore in 

some depth plane wave propagation within a uniform medium ~ i.e. ^efij'^) - ^eff{^) . To that 
end we consider a plane wave solution in the form 
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E(r, co) = E(a)) exp(- / r • k) = E{(o) exp[-i (xk^ + yk^ + z fcj] [ VIII- 3 ] 



which may pictorially represented as 



Therefore 




y 




V co) = V jE(a)) exp(- / r • k)] =- / k • E(r, to) 
V xE(r, co) = V x[e(co) exp(-/ r • k)] = - / k x E(r, co) 



[ Vm- 4a ] 
[ vm- 4b ] 



and the Maxwell's equations formulated in Equation [Vni-2 ] become 
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-/ kxE(f,a)) =-/a)|i.oH(r,co) 
-/ kx H(f ,0)) = / o)e(o)) E(f ,0)) 
-jkE(r,co) = 
-i k • H(f ,(o) = 



[ Vm- 5a ] 
[ VIII- 5b ] 
[ vm- 5c ] 
[ VIII- 5d ] 



Operate through on both sides of Equation [ VIII- 5a ] with the operator "k x " we obtain 

kx[kxE(f,a))]=a)|X(, kxH(f,a)) [VIE- 6a] 

Using the "bac-cab" rule^^ and Equation [ VIII- 5b ] this becomes 

k [k • E(r,co)] - [k • k] E(f ,co) = -w' [i^ 8,^(0)) E(r,co) [ VIII- 6b ] 



or finally 



[k-k] E(f,a))=a)^^io£e£fHE(f,a)) ^ A;^= co^q 8,g(a)) [ VIE- 6c ] 



Substituting these results into Equation [ VIII- 5a ] we obtain 

H(f,(o)= ((0|io)"' k [kxE(f,(o)]= Je,g((o)/|io [kxE(f,co)] [VIII-7] 



so that the wave impedance is given by 

r((xi) = |E(?,a))^|H(?,a))|= Jjlje^^ 



[ vm- 11 ] 



30 



That is a X (bxc)= b(a-c)-c(a- b). 
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Thus, the complete expression for an electromagnetic plane wave propagating in a direction 
k in a uniform medium is given by 



E{r,t) = E((o)exp[-7 (r- k-co r)j 
H(r, r)=[ri(a))]"' [kxE(r,(o)] 



[ VIII- 9a ] 
[ VIII- 9b ] 



ELECTROMAGNETIC INTERFACIAL CONTINUITY CONDITIONS: 



The previous section gives a complete plane wave solution within a particular uniform, 
linear, isotropic medium. The key remaining problem is to find how that solution may be 
extended into a second uniform, linear, isotropic medium. The conditions for extending 
the solution across an interface between two materials are give by consideration of the 
appropriate integral forms of Maxwell's equations — viz. 

d 



|E(r, t) ■ dl = -— II B(r, t) ■ dA 



[ Vm- 10a ] 



^H{r, t) ■ dJ = ^^J{r, t) ■ d A +j-J^D{r,t) ■ d A [VIII- 10b] 

Applying these equations to the small thought loop that spans the interfacial surface, as 
illustrated below 



"Thoughtl'bop 



Medium 1 
[8i'ai' 



i*-^^ Interface 



/ 

AL 



X Medium 2 

[£2 '02' |I2] 



it is seen that Equation [ VIII- 10a ] yields 
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^E{r,t) ■ dl ^{E^{r, t)-E^{r,t)] -x AL =>0 [VIE- 11 ] 

unless B(r, t) is pathologically large over the loop. Similarly, it is seen that Equation [ 
Vlll- 10b ] yields 

^n{r,t)-dl ^{il^{f,t)-H,{r,t)} -T AL =>0 [VIE- 12] 

unless j( r, t) and/or D(r, t) are pathologically large over the loop. 

In words and in general, the tangential component of the electric 
field strength E{r,t) and the magnetic field strength H(r, t) are 

continuous across an interfacial surface between two materials 
unless the electric current density j{r,t), the magnetic flux 
density B{f,t), or the electric flux density D(r, are 
pathologically large near that interfacial surface. 
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The Mathematical Representation of Fields: 

In abstract vector form, the incident field is given by^^ 



. - 1 / ^ ^ [Vin-13a] 



Note: In this figure we have taken the plane of reflection to be identical to the plane of incidence. While 
assumed here for simplicity, this important identity is establish in the analysis below. 

A note on notation: The subscripts _L and 1 1 refer to the polariztion of the electric field taken with respect to the 
plane of incidence. The _L field components are also called transverse electric or TE components and the II 
field components are called transverse magnetic or TM components. 
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the reflected field is given by 

E-f = {Ef - ri 1 k'"'x } exp(-/ k, V.'^'- ?) 



[ Vm- 14a ] 



[ vm- 15a ] 



and the transmitted field is given by 

E'"'"= {e^- k'^^"xHr} exp(-/ k'™?) 

jjtran^ |gtran_^^-l ^ tran ^ j, tra„| gxp (- / it ^ k f ) 

In coordinate form these equations become: 

ginc ^ {Ef y- Ti J- cose,„, x+ sine,„, z]x[Hiry]} exp [-/fc,(-xcose,„, +z sinG^^J] 
H- = {H;ry+ili'[-cose,„, x + sine,„,z]x[Ery]}exp[-/fc^(-^cose,„,+zsine,„J] 

E-f = {Ef y- 111 [cose^^f x+ sine,,, z]x[H-^y]} exp [-ik, {xoos%^^, + z sine J] 
H'^^= {H[ry+Ti-'[cose,,, x + sine,,f z]x[Ef y]} exp[-/fci (xcosS^^.-i-zsine J] 

gtran^ {Er-y- Ti , [-cose,,„x + sin0,,„ z]x[H;r"y]} exp [-ik, (-^cose,,„+z sine J] 
H'-= {H--:y +Tl-2^[- cose,,„ x+ sine,,„ z]x [Ef-y]} exp \-ik, (-xcose„,„+ z sine J] 

Or expanding out the cross-products: 

E"' = {Ery+(TliH7) [cosei„, z + sine^^, x]}exp [-ifcj(-^cosei„, +zsinei„j] 
gmc ^ |H;-y- (n^' E7) [cosei„, z + sine^^, x]} exp [-/ k, (-xcose^^, +z sine^^j] 



[ vm- 13b ] 



[ VIII- 14b ] 



[ vm- 15b ] 



[ vm- 13c ] 
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|Efy + (riiHf) [-cosG^^f z+ sine^^^ x]} exp[-/^i(xcose^^f + zsin0^j] 
H'''= {h;' y - (t|^' Ef ) [- cosG^^f z+ sinG^^j x]} exp[-/ (xcosG^^^ + z sinG^j] 

E«-= {Ery+(Ti, Hir) [cosG,,„z + sinG,„„x]}exp [-ik ,{-xcosQ,^^+zsmQ,j] 
H"'"= {h|^' y- (ti"2 e;;;'") [cosG,,,„ z+ sinG,,,„x]} exp [- j fc^ (-xcosG,,,„+ z sinG, j] 

Applying any kind of continuity conditions at the interface requires that 

G ,= Laxv of Sinus [VIII- 16a] 

K sin0tran= K ^"^^ mo ^-^.w of SttcCC [ VIH- 16b ] 

Applying, in particular, the continuity conditions discussed in the previous section ~ viz. 

[Ex] =[Ei] and [h^] = [h1] [ VIII- 17 ] 

>■ J tang >■ J tang >■ J tang >■ J tang 



at the interface, requires that 

ginc _|_ gref p^tran 

Tl-^cosG,„, [Ef - Ef ] = Ti-^cosG,,„[Er] 

and that 

Hinc I Tjref Tjtran 
II +^\\ 

tIjCosG^, [Hf - Hf ] = Ti,cosG„,„[Hr] 

These two sets of equations yield the Fresnel Reflection Equations ~ viz. 

Ef _ rri-'cosG.„,-ri-'cos6,,J 



[ vm- 18 ] 



[ VIII- 19 ] 



E 



1"' [tI^'cosGj^^+ti^'cosG, 



[ vm- 20a ] 



[ VIII- 14c ] 



[ vm- 15c ] 
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and 



Hr L^iCosei„,+Ti2cose,,,„ 



[ VIII- 21a ] 



Since r] j ' sin = T] 2 sin ^ 



and 



Ef ^ cos0i„, sin0.,,„- cos0,,,„sin0i„, ^ sin(0,,,„-0i„,) 
E"" cos0ij,(, sin0,j3„+cos0(janSin0i„c sin 

inc / 



[ VIII- 20b ] 



h; 



^ cos0ine sin0i„^ - cos0.,,„sin0.,^„ ^ tan(0i„,-0,J 

Hr cos0i„, sin0i„, + cos0,,,„sin0,,,„ tan (0i„, +0, J 



[ Vm- 21b ] 



These equations taken together with first equations from Equations [ VIII- 18 ] and 
[ VIII- 19 ] yield the Fresnel Transmission Equations ~ viz 



Ef 



2cos0i„, sin0^,,„ 



cos0i„, sin0,^^„+cos0,,,„sin0i„^ 



[ VIII- 22 ] 



and 



2 cos0,„,. sinO, 



H|| cos0i„, sin0i„^ + cos0^,,„sin0„^„ 



[ vm- 23 ] 
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0.2 



FAMOUS Fresnel Reflection Curves ( n2/«i= ^i/«2= 7^2/^1 = 1-5 ) 

Hra 
■I 



Gj^j (fraction of -n/S) 






0.2 0.4 



The minimum (zero) in H[|''7h["'' occurs at the Brewster angle where 



0.6 O.S 1 
Gi„ (fraction of tt/S) 



. /^vBrewster, ^Brewste)l 



or 



e Brewster _/t q 
tran - Tt/^-f 



tran y 

Brewste 



inc 



[ Vm- 24a ] 
[ vm- 24b ] 



or (from Snell's equation) 



tane^~= Tii/Ti2 = njn,= Jejei 



[ VIII- 24c ] 
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Total Internal Reflection 

Reconsider Equation [ VIII- 15c ] and use Snell's law to write the exponential factors in 
the form 



E' 



= |E;["y + (ri2H;['") [cose,^^„z+sine^^^„x]}exp ix J^'- ^'sin'0,„^ exp[-/z ^sinGj [ VHI- 25 ] 



When sin9i„^ > ^2/^1 = «2/«i=sin0|^^', E^^ , the solution in medium 2, is 
attenuated! 



incident 
beam 




attenuated 
field 



Reconsideration of Equation [ Vni- 20a ] and [ VIII- 21a ] shows that the magnitude of 
the reflection coefficients are one when sin9:„^ > sin9=!,"' ~ viz. 



Ef [ cose,„,-/Jsin^e,,-(y^)^ ^ 
Ex Lcose,„, + /Jsin'e,„,-(v^) J 



= exp-j-j 2 tan 



r 

-i| V 



sin'ei„,-(fc2/^i) 



cos9=„ 



[ Vm- 26a ] 



and 
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H 



ref 



h: 



= exp 



- 7 2 tan ^ 



COS0, 



[ Vm- 26b ] 



PARALLEL PLATE WAVEGUIDE ; 

Consider the propagation of a plane wave between two parallel perfectly conducting planes. 





Perspective view 



Side view 



Combining Equations [ VIII- 13c ] and [ VIII- 14c ], the electric field strength of the TE 
wave in the region between the plates may be written 

Ex = y E"' exp(ix k, cos0i„,) + Ef Qxp(-ix k, cos0,„,)] exp(-/ z k, sinOi^J [ viII-27 



At X = the field parallel to the surface of a perfect conductor must be zero so that 



E^^'^=-E"'' and, therefore, 



Ei = yE"'' exp^/x cos9i„^j- exp^-/xA:iCOS0,„^j exp(-/z ^ sinGj^J 
= y 2i E"" sin ^xk^ cosGj^^. j exp(- i z P) 



[ vni-28 ] 
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where P= sinGj^^ . At the upper surface ~ /. e. x = d — the field parallel to the surface of 
a perfect conductor must also be zero so that 

dk^cose^^^=nn where n= 1,2,3,... [ VIII-29 ] 

and, therefore, 

Pte„ = K sinGi^, = Jk^- kl coiQ-^^ = - (n njdf wheren= 1,2,3,... [Vni-30] 

which is the dispersion relationship for TE waves in a parallel plate 
waveguide with "cutoff" frequencies at 

tof = {n Ti/d) (Ji^Fj) wheren= 1,2,3,... [ vni-31 ] 
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Again combining Equations [ VIII- 13c ] and [ VIII- 14c ], the electric field strength of the 
TM wave in the region between the plates may be written 
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E||= z(rii cos0i„,) {h7 exp(ixk^cosQ,^^)-ll\f oxp(-i xk^cosd^^^)^ex^-i zk^sinQ,^^) 

+x(TiiSinei„J {h"' exp(/x k, cosGJ + H,'' exp(-/x k, cosGi^J} exp(-/z k, sinG^^J ^^^^^ 



At X = the field parallel to the surface of a perfect conductor must be zero so that 

Iref Tjini 
II - ^1 



H[|^^= H"'' and, therefore. 



E„= z 2 i Hir (ill cosGi„,) sin(x k, cosGi„,) exp(-iz p) 

/ ^ / \ / [VIII-33] 

+x 2 H- (t1i sinG,„,) cos(x k, cosG,„,) exp(-/z P) 

where P= k^ sinGi^^.. At the upper surface ~ i.e. x=d — again the field parallel to the 
surface of a perfect conductor must also be zero so that 

dk^cosQ-^^=nn where n= 0,1,2,3,... [ VIII- 34 ] 

and, therefore, 

= k^ sinGi„, = yjkf-kfcos^e^^ = ^k\-{n njdf wheren= 0,1,2,3,... [ VIII-35 ] 

which is the dispersion relationship for TM waves in a parallel plate 
waveguide. 
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Note that the TMq mode is a bona fide mode of propagation which does not 
have a "cutoff" frequency! 

Dielectric Slab Waveguide ; 

Consider the propagation of waves "trap in" or "guided by" a dielectric slab of thickness d. 

X 

A 
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In its full generality this is moderately complicated problem, but a rather simple ray optics 
model of the propagation is sufficient to yield dispersion relationships for the various 
possible modes of propagation. To obtain such relationships, consider the total internal 
reflection of a sequence of plane waves as illustrated below. 





Ginc ; 

y 


f k 








1 


\ ^^^^^^ 






X d 








^, ~^ 


<^ < 

s 




z 





T 



In order for the multiplely reflected wave to be self-consistence the following, relatively 
obvious, phase condition must hold:^^ 

A(^^^^ +A(^^^o+ 2 Ai J cos0i„^ =m27C where m = 0,1,2, 3,... [ VIH- 36 ] 

where A(|) and A(|) are, respectively, the phase shifts associated with the reflections 
at the upper and lower dielectric boundaries. 

For TE-modes of propagation Equation [ VIII- 26a ] gives the phase shift at the 
boundary (called in the trade the TE Goos-Hanchen shift) and Equation [ VIII- 36 ] 
becomes 



This equation is a direct generalization of Equations [ V-29 ] and [ V-34 ] which figured in our analysis of 
parallel plane waveguides. 
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cos 9:, 



cos 8:. 



or 



— - — Jtif-n =m7r+tan ' — ' 

V ^0 y 



+ tan 



[ VIII- 37a ] 
[ VIII- 37b ] 



where sinGj^^^ = ^fkj = njn^ (n is the effective index of the propagation mode). For the 
symmetric case (i.e., = nj, the self-consistence relationship for the TE modes is given by 



cos 0. 



tan 



^0 J cosGj^j 7t 

m — 

2 2 



[ VIII- 38 ] 



where = 0)/c= Ih/Xq- This is a transcendental equation in the single variable cosG^^^. 
Its solutions yield the allowed bounce angles, (Gj^j,)^^ , of possible modes and, hence, the 
allowed propagation constants since P= sinGj^^ The left and right sides of this equation 
may be plot as a function of cosGj^^, with kQd=n^ 2% [clfko) and sinG°"= nju^ as 

a parameters. The intersections of such curves yield the allowed bounce angles as 
illustrated below 
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LHS and RHS of Equation [ VIII- 38 ] 
for n, {d/\)= 0.5 and cos0i:'= 0.5 




0.1 0.2 0.3 0.4 0. 



LHS and RHS of Equation [ VIII- 38 ] 
for n, {dl\)= 1.5 and co%%Zl= 0-5 
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'inc 



LHS and RHS of Equation [ VIII- 38 ] 
for n, {d/\)= 1.0 and 0080;:'= 0.5 



0.5 




LHS and RHS of Equation [ VIII- 38 ] 
for n, {d/\)= 2.0 and 0080;;'= 0.5 



0.5 
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1974, Kogelnik and Ramaswamy ^4 developed a convenient formalism for treating slab- 
waveguide problems. First they introduced three new waveguide parameter ~ viz. 

The normalized frequency/slab thickness parameter V= k^d •^n\ - n\ [ VIII-39a ] 

The normalized waveguide index parameter b={n^ - n\^l{n]- n^^ [ VIII-39b ] 

The normalized waveguide asymmetry parameter «-(«s~"c)/("f~"s) t Vni-39c ] 

They then showed that Equation [ Vin-37b ] could be written 

V JT^ = m 7U+ tan"' ^bl{\-b) + tan"' ^{a + b)/{l-b) [ Vni-39c ] 

which can be used to generate the following family of curves: 



H. Kogelnik and V. Ramaswamy, Appl. Opt. 13, 1857 (1974). 
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IX. Optical Pulse Propagation 

The Electromagnetic Nonlinear Schrqdinger Equation; 

We begin our discussion of optical pulse propagation^^ with a derivation of the nonlinear 
Schrodinger (NLS) equation. To that end, we recall Equations [ VII-23 ] and [ VII-23 ] 
from the early lecture set entitled Nonlinear Optics I — i.e. 

V' E(r,co)- V(V •E(r,a)))+(to7c'eo)l(a)) •E(r,a)) = -^lo co' P*'''''(r,a)) [ IX-1 ] 
V (g(a)) • E(?,a))] =- V •P^'"^'(?,a)) [ IX-2 ] 

In this treatment we will confine our attention to wave propagation in uniform, isotropic 
optical materials ~ viz., glass fibers. For such materials, we can write 

p(^)(f,a))=8^(f,a))E(f,a)) [IX-3] 

3 - 2 

where 8f^(r,0)) = —eg xl^L|E(f >0>)| and, thus. Equation [ IX-I ] simplifies to 

V' E(f ,©) +(/:o7eo) [eH +eNL(f ,co)] E(f ,©) = [ IX-4 ] 

where kf^=o)/c.^^ 

To proceed, postulate that this nonlinear Helmholtz equation can be treated by separation 
of variables methods. In particular, we are looking for a time-localized solution (a 



An excellent reference on this subject is Govind P. Agrawal's Nonlinear Fiber Optics, Academic Press (1989) 
ISBN 0-12-045140-9. 

In this simplification, we have taken V •£(?, Co) = -[e(o)) -8 NL(f , »)] E(r,a)) • V eNL(f ,») ~ . 
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pulse) with a relatively narrow frequency spectrum (or "group" of frequencies centered on 
a frequency co^.^. Thus, we assume a seoaration of variables solution 

E(f ,co)= F{x,y) G{z,(0-(oJ exp( - / [ 1^-5 ] 

where [3^^ is a wave or propagation number to be associated with to^.^ and, thus. Equation 
[ IX-4 ] becomes 

where k= J[8(o))+8^(r,to)]/eo . In the linear problem would be the "separation 

constant," but in this case we will need a bit more elaboration. Nevertheless, we shall 
assume that we can find a set of functions F(x, y) and values that satisfy the equation 

+ + -P ]f=0 [IX-7a] 

sothat p_,-2p,,|^ + [I5^-py G=0 [IX-7b] 

To use perturbation theory, we first reduce Equation [ IX-7a ] to a solvable linear problem 
by writing 

= kl [£((o)+8NL(f ,co)]/8o= [n(a))+A n{(ii)\ 
~ [n^(co) + 2 n (co) A n (co)] k^ 

P' = [P+AP]'»P' + 2PAP [IX-8b] 



[ IX-8a ] 
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eNL(r,«) 3 xi!L(«) If;/- xp 
where A n (coj = — = ^--^ |Ji,(r,(Oj| . Thus, to first order we need to solve 



2 n (to) 8 n (to) 



the linear equation 



a^F r 2 .2 r21p . 



[IX-9] 



which taken together with appropriate boundary conditions defines the linear eigenvalue 
problem for propagation in the medium where the functions F are the eigenfunctions and 
the values |3 are the eigenfunctions. In the previous lecture set — i.e., VIII. Guided Waves 
in Planar Structures — we found a set of eigenfunctions and eigenvalues appropriate to the 
dielectric-slab guidewave propagation 



Following an earlier discussion, we now presume that 



oz 



i.e., we take 



theslowly vary amplitude or paraxial approximation ~ so that Equation [ IX-7b ] 
reduces to 



ac 



where J5'-p'^»p'-p'^+2p Ap and, thus. 



ac 



[IX-lOa] 



-'■2P-3- + [P'-PcjG+2pApG=0 



[ IX-lOb ] 



For equation [ IX-7a ] to be completely satisfied in first order, we must have 
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_ nk^ JJ An |F {x,y)f dx dy 
P jj |F {x,y)( dx dy 



8P 



(3)1""/" '\| jj| v '-'/| "'^"y 



Since we are assuming that propagating pulse has a relatively narrow frequency spectrum, 
it is reasonable to use a Taylor expansion around co^.^ for p (co) ~ viz. 



PH==Pctr+(«-cojPi+^j((o-(oj'P2 + ^j(co-(aj'p3 + , 



[IX-12] 



where - 



d%c^) 



day' 



Near the dispersion minimum in glass fibers {i.e. X~ 1.3-1.6|i.m ) we may, to very 
good approximation, stop with the quadratic term and write 

r 1 2 1 

P^-|3L»2|3,,[j3-|3j = 2p,,^(a)-Q)jp, + -(Q)-a)J p,j . [lX-13] 



In this approximation. Equation [ IX- 10 ] becomes 

1 



1^ + / [(co-to J p,+i(a)-a) J'p.j G+i^ Ap G = 



[IX- 14a] 



If we take 



P A o ^0 (3) J 



]\F{x,y)( dxdy 



J |F(x,3;)| dxdy 



[ IX-15 ] 
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Equation [ IX- 10 ] in the frequency domain reduces to 
dG 



+ 1 



1 2 

(a)-o)Jpi+-(to-toJ % 



g=-/y|g|'g 



[ IX- 14b ] 



which implies the following time domain equation for the pulse envelope: 
BG{z,t) _ BG{z,t) / _ d'G{z,t) a - 



dt 



■ + ^G{z,t) = -iy\G{z,tfG{z,t) [IX-16] 



The last term on the left hand side has been added to incorporate the effects of various 
possible loss mechanisms. 

Next we transform into a coordinate system which moves with the "group" ~ what might 
be called the surfer's coordinates of the pulse ~ i.e., {t,z}^{i,q} where 
T= Pi (z- Zq) and q= z- Zq= z-^i^ t^. In term of these surfer's coordinates. 

Equation [ IX-16 ] becomes^'' 



dq 2 



= - ^ G(G, t) + ^P, - i Y|G(G,x)f G(^, t) 



.2 



[ IX-17 ] 



Obviously, if we omit all of the terms on the right hand side of this equation so that 
dG{q,x) 

— — = 0, we would have the ideal situation wherein a pulse of any shape 
propagates forever without changing shape at a velocity v = Pi~' —the group 



Since 



|c^,,,).i5M|.lGfe.)| = lGfe.)W.AG,,.)K) 

iG(...) = |GM|.^Gfe.)| = iG(;..)(0).lG(;..)(l)- 
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velocity. In treating the less-than-ideal situation, we will initially neglect the effects of the 
loss (first) term and confine our attention to the competing effects of the dispersion 
(second) and nonlinear (third) term . We cast the lossless version of Equation [ IX- 16 ] 
into a normalized, standard form by introducing 

JPq Ld Ax Ax 

where Ax is the width of the pulse and Pq is its peak power. 

Thus, we, at last, obtain the standard form of the nonlinear Schrodinger (NLS) 
equation 

du(Z,T) . . d^u(Z,T) , ,,2 . , 

/— ^-sgn(|3,) ^ + N'\u{Z,T)\ u{Z,T)= [ IX-18 ] 

where N' = LJL^^= y P, Ax' 

PULSE Solutions of linear Schrodinger equation 

If in Equation [ IX- 10 ] we set G[z,t) ^ U[z,t) and neglect the loss and nonlinear terms, 
we see that U(z, t) satisfies the following differential equation: 

— lXz,t) + v-;-U{z,t)=ib ^iXzj) [IX-19] 



For minimal dispersion ~ viz. n b = — — ■ — — 



2! da 

becomes 



»0 -- Equation [ IX-19 ] 



|- U(z, t) + v"' I- U(z, t) = [ IX-20 ] 

dz at 
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which is the basic wave equation for minimally dispersive media with the general 
solution 



where v„ = 



^^P(co) 






10=0) Q _ 



U{z,t) = u(z-v^t) 
is the group velocity of the pulse. 



[ IX-21 ] 



When limitation to first order dispersion is an adequate approximation, Equation 
[ IX- 19 ] expressed in surfer's coordinates becomes 



[ lX-22 ] 



where b= ^J^- 



Amazingly, this pulse dispersion equation is the, so called, parabolic 
equation that we saw earlier in connection with beam propagation -- viz. 
the paraxial wave propagation equation. 

Solution of Pulse Dispersion Equation 

For convenience, we restate here Equation [ IX-22 } the first-order pulse dispersion 
equation ~ viz. 

b£-,\J{q,x) + i^^U{q,x) = 0. 



Let us write a Fourier transform for this modulation in terms of "surfer time" ~ i.e. 
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iXq, x) = J U((;,(o) exp( j(o x) dco [ IX-23a ] 

where 

IJ{(;, «) = ^ J U(<;> t;) exp(- / CO x) Jx . [ IX-23b ] 



Thus, the pulse dispersion equation reduces to an ordinary differential equation - viz. 

i Y U(<;'(o) = ^ «^ U((;, 0)) [ IX-24 ] 

for the Fourier transform and thus we have the simple solution 

U((;, co) = U( 0,0)) exi(- ib03\) . [ IX-25 ] 

Thus, we see that the dispersion changes the phase of each spectral component 
of the pulse by an amount that depends on the frequency and the propagated 
distance. The general solution may be written as 

x) = J U(0,a)) exij/ (co x- b to' c;)] Jto [ IX-26a ] 



where 



U(0,to)= — Ju(0,?- ?o) exp[-/ to(r- ^o)] d{t- t^) ■ [ IX-26b ] 



27U _ 

Let us suppose that we have a "chirped" Gaussian at c; = ~ i.e. 
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U(0,r- exp 



r 



l-iC 



t-tA 



V '^0 y 



[ IX-27a ] 



so that 



U(0,a)): 



Xn CO I 



27l(l-/C) ^'^^ 2(1- jC) J 



[ IX-27b ] 



where C > characterizes an "up-chirp" and C < a "down-chirp" pulse. 



A voltage 




A Gaussian pulse with a frequency "down-chirp" 



Inverting the transform, we see that 
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, , r 

u((;,t) = 

or rationalizing 

r 



exp{ 



{l-iC)x' 



1 



xl + i2bq{l-iC)\ ^ 2[xo + /2Z7(;(l-iC)]j 



1 



2x; 



2 / 



l+- 



+ 



^0 J 



1+ 



^0 y 



xexp ■ 



i 




1 + C') 




2x1 


'( 


''0 / 


2 


I '^0 





Hence, the pulse width broadening factor at a given position z is given by 



1 + 



+ 



2bq 

V ^0 y 



2 ^ .2 



V 



+ 



where Ljj = xl/2b. 



[ IX-28a ] 



[ IX-28b ] 



[ IX-29] 
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down-chirp" 



C = 



"sero-chirp" 



0.5 1 1.5 2 

RELATIVE DISTANCE (z/Lp) 



The spatial evolution of the pulse width of chirped Gaussian pulse 



(For "normal dispersion" -- i.e. 



>0) 



Pulse Solutions of "Negligible Dispersion" Nonlinear Schrodinger Equation 



If group velocity dispersion can be neglected. Equation [ IX- 17 ] reduces to 



2 - 



M=_|Gfcx)-,y|Gfex)-G(,,x) 



[ IX-30 ] 



and if we take (^q, x)^U{q,x) exp (-a<;/2) we obtain 
dU{q,x) 



=- i Lnl exp(-a(;) | U{q,x) \ U {q,x) 



[IX-31 ] 
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where the characteristic nonlinear length is given by L^^ = (y ^) - A solution to this 
equation is readily obtain in the form 

[ IX-32a ] 



where 



and 



U{(;,x) = U{0,x)^xp[i^^{q,x)] 

^nl{<;,'^) = -[l- exp(-a(;)] [aL^L]"' exp(-a(;)| U (0,t) 
= -[<icJ^m]\U{0,x)\" 

Geff=[l-exp(-a(;)] a"' 



[ IV-32b ] 
[ IX-32C ] 



This interesting result show that so called "self -phase modulation" or SPM gives rise to an 
intensity dependent phase shifted or chirped pulse which remains constant 
in shape as it propagates. The instantaneous optical frequency shift is given by 



i^) = = -['ieff/LNLj 



ax 



[ IV-33 ] 



Note that the pulse spectrum is "red-shifted" on the leading edge of a pulse and "blue- 
shifted" on the trailing edge of the pulse. If we suppose the initial pulse to be a super- 
Gaussian of mth-order ~ i.e., 



C;.(0,T) = exp 



r i+i 



iC 



[ IV-34 ] 



~ then the instantaneous optical frequency shift would be given by 



5a) (t) = 



Seff 

v^NL y 



/ N2m-1 



vAXy 



exp 



vAXy 



[ IV-35 ] 
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A Brief Discussion of Solitons 

We return briefly to Equation [ IX- 18 ] — the standard form of the nonlinear 
Schrodinger (NLS) equation - i.e., 

du(Z,T) . , d^u(Z,T) , ,,2 . , 

i — ^ - sgn(|3 J ' + N' \u{Z,T)\ u{Z, T)= [ IX-18 ] 
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where M(Z,r) = -^^, z=^=<;^-pj, T=J~2^ and N'= LJL^= yP,Ax'l\^, 
^Pq Ld At At 



In soliton analysis the most common fomi of NLS equation is the following: 
dU{Z,T) d^U{Z,T) 



dZ dT' 

'y At^ 



+ \U{Z,T)\U{Z,T) = [IX-36] 



where U{Z,T) = N u{Z,T)^G{(;,x) 



P2 



If A'^ > 1 the following fundamental soliton will propagate undistorted for an arbitrary 
distance: 

U{Z,T) = sech{T) exp(j Z/2) [ IX-37 ] 

If N> 2 the following second-order soliton will propagate undistorted for an 
arbitrary distance: 

4 [cosh(3 T) + 3exp (/ 4 Z) cosh(r)l 

U(Z,T) = r , \ , , 7^ exp(i Z/2 [ IX-37 ] 

^ ^ [cosh(4r) + 4cosh(2r) + 3cos(4Z)] ' ^ 

2 16[cosh2(3r) + 9cosh'(r) + 6cos(4 Z) cosh(3r) cosh(r)] 

r i^' ^)\ ~ r T2 [ IX-37 ] 

[cosh(4 r) + 4 cosh(2 T) + 3cos(4 Z)] 
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Two SoLiTON Collision illustrations 



The following are time-lapsed illustration of the propagation and collision of two 
solitons ~ plot here is 

. . ,^ 3 + 4cosh(2x-8?) + cosh(4x-64?) 

U[X, t) = 12 — — ;^ ^li 



[3cosh(x-28?) +cosh(3x -36?)] 



8 



time = -1 .00 




-20 



-10 



8 



time = -1 .00 



li 




position -20 



-10 



position 



time = -0.25 


8 






6 






4 






1/ 





-10 



position 
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